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Abstract. We study formally and rigorously the bifurcation to steady and time-periodic 
states in a model for a thin superconducting wire in the presence of an imposed current. 
Exploiting the PT-symmetry of the equations at both the linearized and nonlinear levels, and 
taking advantage of the collision of real eigenvalues leading to complex spectrum, we obtain 
explicit asymptotic formulas for the stationary solutions, for the amplitude and period of the 
bifurcating periodic solutions and for the location of their zeros or "phase slip centers" as 
they are known in the physics literature. In so doing, we construct a center manifold for the 
flow and give a complete description of the associated finite-dimensional dynamics. 

1. Introduction 

One of the natural applications of superconducting is to exploit their infinite conductivity 
to transmit electric currents. The goal of this paper is to analyze a number of asymptotic 
problems that arise in the study of such current transmission through a wire. We consider a 
simple canonical problem, in which the superconducting portion of the wire is of a finite extent 
—L < X < L. It is assumed that a normal current / is fed into the the wire at its left end. It is 
known that under the right circumstances, for example for a temperature T that is sufficiently 
small, the current in the wire itself will be in part normal and in part superconducting. This 
coexistence of two types of currents in the wire is called a resistive state. 

The resistive state in superconducting wires received some attention by physicists who ob- 
served a variety of phenomena that are unique to this situation. To review these observations, 
which involve oscillatory, that is to say, inherently time-varying behavior, it is standard prac- 
tice to use the time-dependent Ginzburg-Landau model (TDGL). For a three-dimensional wire 
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occupying say a thin cylindrical region D with axis of length 2L centered on the x-axis, this 
system can be written in non-dimensional form as 



(cf. [T7]) where ip : D ^ Cis the Ginzburg-Landau parameter whose square modulus measures 
the density of superconducting electrons, A : M'^ ^ M'^ is the magnetic potential whose curl 
measures the effective magnetic field and : — M is the scalar electric potential whose 
gradient represents the electric field. The parameter F is proportional to Tc — T, where T is 
the actual temperature, and Tc is the phase transition temperature in the absence of external 
currents (i.e. in the case / = 0). In (11. 2p . * denotes complex conjugation and the right- 
hand side represents the sum of normal current (with associated ohmic conductivity a) and 
supercurrent. We note that the TDGL is invariant under the gauge transformation 

{ip, 0, A) — i> ('06*^, (p — gt,A + Vg) for any smooth scalar function g = g{x, y, z, t). 

To pursue an appropriate three-dimensional analysis of the problem of forcing an applied 
current into a wire, one would then have to impose inhomogeneous boundary conditions on the 
normal component of the normal current and couple the system above to a Maxwell system 
on the exterior of D. This is not the direction we will follow; instead we adopt the model 
favored for many years in the physics literature on the subject, e.g. fS], and view the wire as 
a one-dimensional object. Before leaving the higher dimensional setting we should comment, 
however, that an interesting two-dimensional study of the stability of the normal state with 
an applied electric current can be found in [1]. 

In such a one-dimensional model, the exterior problem is typically ignored as a lower- 
order effect and so all unknowns are taken simply to be defined along the wire as functions 
of X and t only. Furthermore, through the gauge choice gx = —A, one can eliminate the 
magnetic potential A completely from the system. Then using the remaining freedom in the 
i-dependence of g, one can insist on the convenient normalization 



(1.2) 




iP^ + i(f)ilj = (V - iAfip + (F - for (x, y, z) G D, 

V X V X A = -a{At + V0) + i (^^* - ^^^*) - A for (x, y,z) eD 




0(0, t) =0 for all t > 0. 
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Under these assumptions and gauge choice, (II .ip reduces to 

gll I (1.4) + = V'xx + - 1^1^^ for - L < X < L, t > 0. 

Regarding the fate of the second equation (ll.2p . note that necessarily the divergence of the 
total current, that is, the right-hand side, must vanish. In one dimension, this condition 
however implies simply that the total current is a constant. Therefore, since we are specifying 
that the current at the endpoints of the wire is purely normal and equal to J, we arrive at the 
relation 

gis] (1.5) ^-{ijiljl-ilj^ij*)-a(j)^ = I for -L < X < L, t > 0. 

Since the natural setup is for the wire to be connected at its endpoints to a metal exhibit- 
ing normal conductivity, we supplement the system f ll.3p - fll.5p with homogeneous Dirichlet 
boundary conditions on the order parameter 

gl5] (1-6) V^(±L,t) = 0, 

along with initial conditions on ip. The general nature of our results apply also to other 
homogeneous boundary conditions, including in particular the case of homogeneous Neumann 
boundary conditions. 

As long as the temperature is sufficiently high, that is F is low, the wire is in the normal 
state. This means that ^ = 0, and the current in the wire is purely ohmic, i.e. 

glT] (1.7) <p. = -I/cr. 

As the temperature T is lowered, one reaches a critical value, determined by a curve F = Fi(/) 
where the normal state loses its stability and a nontrivial superconducting state might emerge. 
Two such states were discovered in the early 1980 's. The first one is a stable stationary 
solution, reported by Langer and Baratoff [9]. The notion of 'stationarity' requires some care 
here. Let us express the order parameter in polar form ip = /e*^; then the gauge invariant 
quantities f{x,t), q{x,t) := Xx{x,t) and 6{x,t) := Xt{x,t) — 4>{x,t) converge to stationary 
functions fo{x), qo{x), 9o{x). On the other hand, Langer and Ambegaokar ([IT] (see also Ivlev 
and Kopnin [8j) found in their numerical simulations that for some values of the parameters 
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Figure 1. 

(/, r), the normal state bifurcates into a state where the order parameter oscillates in time. 
This periodic behavior is particularly interesting in light of the dissipative nature of the TDGL 
model (ll.4p . but it is made possible by the presence of the applied current which effectively 
disrupts the gradient-flow structure of the system. 

For the duration of this investigation, we fix L = 1 so that the wire occupies the interval 
[—1,1]. We also set the conductivity a = 1 in order to focus on the different kinds of states 
that emerge at different points in the (/, F) plane. 

In a recent work [13] the authors used numerical simulations and some analytical arguments 
to identify a more elaborate phase transition picture. In particular the curve Fi(/) was shown 
to be associated with an interesting spectral problem. It was also shown that there exist two 
critical currents and Jc that play an important role in the system behavior. The curve Fi(/) 
is depicted in Figure [H Specifically, the normal state (N) is only stable for F < Fi(/). When 
the temperature is decreased and F increases past Fi(/), the (N) state become unstable. If 
/ < Jfc ~ 10.92, then the (N) state bifurcates into a stationary (S) state. On the other hand, 
the bifurcation branch to a stationary state is unstable for Ik < I < Ic, while if /> Jc ~ 12.31, 
the (N) state bifurcates to a stable time-periodic (P) state. 

To understand the phase transitions described above, and in general the solution to equa- 
tions fll.4l) - fll.5p . it is useful to study a number of mathematical problems: 
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(1) What is the meaning of the special current value Ic where the bifurcating state switches 
from a stationary one to a periodic one? 

(2) From the viewpoint of self-adjoint operators arising in the time-independent case, 
the answer to the first question involves an unusual spectral problem involving the 
collision of two eigenvalues. Moreover, near the critical value 1^ the spectrum of the 
underlying operator changes its nature, from real to complex. Therefore it is desirable 
to understand this spectrum near the special value Ic- 

(3) What is the nature of the bifurcating branch near the curve ri(/)? In the language of 
dynamical systems, we ask what is the geometry of the center manifold there? This 
question, in fact, involves a number of issues. For instance, is the bifurcation branch 
stable and what is its shape? In addition, we point out that in the periodic case there 
are points in space-time where the order parameter if) vanishes. Such points are called 
phase slip centers (PSC's) since the phase 'exploits' the vanishing of the amplitude to 
have a discontinuity there, thus relaxing large accumulated phase. Thus one of the 
relevant questions would be to identify these points. 

The loss of stability of the normal state is studied through the linearization of equations 
fll.3p - fll.6l) around the normal state 

"jlir] (1.8) ^" = 0, (t)=-Ix. 

It is convenient to express the solution ip of the linearized equation in the form ip{x,t) = 
u{x)e^^~'^^^; then we obtain that u{x) is the solution of the following spectral problem 

gll3| (1.9) Mu := u^^ + ixlu = -\u, m(±1) = 0. 

Clearly the stability of the normal state is determined by whether F is larger or smaller than 
ReAi where Ai refers to the eigenvalue of the operator M having smallest real part. In the 
next section we shall therefore consider this eigenvalue and examine some of its properties for 
small I and for large /. The results of this section will help us in identifying the critical value 
Ic- Then, in Section 3, we will examine in more detail the leading eigenvalue of M for I values 
near I^. In Section 4 we construct stationary solutions to fll.4p - (11.61) for I < I^. using formal 
asymptotic expansions and multiple time-scales. In Section [5l we construct periodic solutions 
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using the same methods in the regime I > Ic where the spectrum of M has become complex. 
We then make these calculations rigorous in Section 6 by constructing the center manifold for 
the solution immediately after bifurcation and studying the O.D.E.'s which govern the flow 
on the center manifold. 

A novel aspect of the rigorous analysis in Section 6 is that both the linearized and the full 
nonlinear system admit what is called PT-symmetry; namely an invariance under the joint 
transformation of x ^ —x and complex conjugation. This type of symmetry has been the 
focus of a number of recent investigations (see e.g. [Zl HI El [12] ) . In particular, in the analysis 
of periodic bifurcation we make strong use of this symmetry to reduce the dimension of what 
turns out to be a four- dimensional phase space to a planar system exhibiting a standard Hopf 
bifurcation. Furthermore, we are able to then go further and describe bifurcation for the 
full four-dimensional system involving possibly non-PT-symmetric solutions, proving that, at 
least in the vicinity of the normal solution ^ = and for values of the bifurcation parameter 
that we study, solutions generically converge, up to fixed complex rotation, to the manifold 
of solutions exhibiting PT-symmetry, and thereafter to the stable, PT-symmetric periodic 
solutions arising through Hopf bifurcation within that manifold. However (Theorem 16. lip , we 
also exhibit through direct calculation the existence of unstable, non-PT-symmetric periodic 
solutions in the same vicinity, i.e., persistent solutions that do not converge to the manifold 
of PT-symmetry. This shows that the observed generic convergence to PT-symmetry is the 
result of detailed local dynamics on the center manifold about the normal state and not, say, a 
global principle associated with a decreasing Lyapunov functional. In particular, there might 
exist attracting steady or periodic states far from the normal state that are not PT-symmetric, 
an intriguing possibility to keep in mind in further investigations. 

Finally, in Section 7, we show that our rigorous construction of periodic solutions leads to 
a proof of the appearance of the phase slip centers, that is, periodically appearing zeros of the 
order parameter. 



resistive state in a superconducting wire 7 

2. The spectrum of the canonical PT-symmetric problem for small / and 

for large i 

In this section we examine the spectrum of M defined in equation (11 .Qp . The operator 
M is not self-adjoint, of course. On the other hand it enjoys a symmetry that is called PT. 
The letter P stand for parity, i.e. transforming x —>■ —x, while the letter T stands for time 
reversal, i.e. complex conjugacy. One readily checks that under this pair of operations, M is 
unchanged. 

PT-symmetric spectral problems seem to have been little studied until quite recently. In 
one of the earliest works on this subject (in a physics context). Bender and Boettcher [2] 
considered a canonical PT-symmetric operator on the entire real line and observed through 
numerical simulations that the spectrum is real. In the case of a finite interval as in equation 
( 11. 9p above, the situation is more involved. In particular, we will formally demonstrate the 
appearance of complex eigenvalues for / large with corresponding eigenfunctions possessing 
an internal layer. 

We look first in the case of small /. When / = 0, the spectrum is of course real, and can, 
in fact, be written down explicitly 

(2.1) Afc(J = 0) = (7rA;/2)% fc = l,2,... 

The PT-symmetry of M implies that if {X,u{x)) is a spectral pair of eigenvalue and eigen- 
function, then {X* ,u*{—x)) is also a spectral pair. Langer and Tretter [12j have shown that, 
as long as the eigenvalues of a PT-symmetric problem are simple, the spectrum is a smooth 
function of the problem's parameters. In our case it implies that since the eigenvalues are well- 
separated at / = 0, they are smooth functions of / at least for I small. However, this implies 
that the eigenvalues must remain real also for / positive (but small), since a real eigenvalue 
can become complex only by splitting into a pair of eigenvalues (by the PT-symmetry). 

What happens when / is increased? As long as the eigenvalues do not collide, they remain 
real. We now show formally that in fact the eigenvalues must collide for some I and establish 
an asymptotic formula for the leading (complex) eigenvalue when / is large. Precise and 
rigorous asymptotics for this spectral problem were carried out by Shkalikov [T3]. His work 
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(see also [12]) was performed in the context of the Orr-Somerfeld equation in hydrodynamic 
stabihty theory and makes use of a change of variables leading to an Airy-type equation. We 
present the following formal calculation with the hope that it makes the emergence of the 
complex spectrum more understandable. Note also that the formal method given here applies 
in more general circumstances, whereas the exact solution of [14] is specific to the exact form 
of the equations under consideration. 

To study the spectrum as / becomes large, it is convenient to introduce a small positive 
parameter e, and then write I = e^^. It is clear that the eigenvalues must also be large to 
balance the large 'potential' ie~'^x. We therefore write to leading order 



The eigenvalue problem (11.91) can be written to leading order as u^x + s^'^Qu = 0, where Q : = 
x+jS—ia. Crudely setting e = 0, we obtain a formal eigenvalue problem {ix — {a+il3))tp = for 
the multiplication operator ix, which evidently has only pure imaginary, essential spectrum. 
From this we may conjecture that the spectrum of L becomes complex as / ^ oo (e: — 0); 
however, this limit is very singular and must be examined in more detail (indeed, on the whole 
line, the spectrum is real for / large [8].) 

Consider first the case where the potential Q does not vanish for x in the interval [—1,1]. 
For example, this would occur if a 7^ 0. If Q 7^ then there is no turning point in a standard 
JBKW expansion of equation (11.91) . Therefore, we seek in this case an asymptotic expansion 
of the form 



^ (2.2) 



^ (2.3) 



u{x) = e' 



.iS{x)/e 



Substituting (12. 3p and ( 12. 2p into (11. 9p we get to leading order 



^ (2.4) 



We use the term JBKW instead of WKB since this expansion method was introduced by Jeffries in 
1923, three years before it was rediscovered by Wentzel, Kramers and BriUouin, who are also now ordered 
alphabetically. 
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Integrating the last equation from —1 to a; gives 

(2.5) = ±1^1/2 ((a; + ^ _ ,«)3/2 _ (_i + ^ _ ,^)3/2^ 

For future reference we introduce the notation 
^ (2.6) S+{1) = ReS + ilmS. 

The general solution to the equation (11 .Op is (to leading order) 



^ (2.7) m(x) ~ Ae'^^(^)/" + Be'^ ^^^/^ 



ml? 



ml 



Substituting this solution into the boundary conditions at ±1, and seeking a pair of nontrivial 
coefficients A, B leads to the following complex-valued equation 

(2.8) ^g.5-(i)/._ 

In particular the magnitudes of the two sides of equation (12.8p must be the same. Equating 
the absolute values, and using the notation (12. 6p gives 

ImS/e _ -ImS/e 



19] (2.9) e'"^^/" = e 



Therefore, a regular JBKW expansion without a turning point is feasible only if Im S" = 0. 
Recalling the definition of Im S, the last condition on it implies (for some real number x) 



^ (2.10) {l + (3-iaf/^-{-l+(3-iaf/^ = Xi~^/^- 

We now show that equation (12.101) holds for any a if /3 = 0. Setting /3 = 0, we write 
^ (2.11) - I - ia = pe'^-^^''\ I - ia = pe'^-^'\ 



Here p = a/T+c?. Substituting this into (12.10P with P = 0, and defining 7 = 3p/2 gives 
i43] (2.12) p^^H~^/^G = xi^^\ G = V2(cos7 + sin7) . 

This proves the assertion above with 



i45] (2.13) X = V2(l + a^)^^^(cos7 + sin7) 
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It remains to show that there exist values of a for which equation fl2.8p is solvable. Since 
equation fl2.8p holds if and only if ImS* = and sin (Re 5") = 0, we obtain the condition 

s47| (2.14) Re5 = en7r, 

where n is an integer. The calculation above for x together with equation fl2.14|] imply 

3/2 

s49 I (2.15) (l + a^) (cos7 + sin7) = ^nvr. 

3 

Consider the left-hand side as a function of a (recall that 7 also depends on a through the 
relation (12. lip ). When a tends to 00 then 7 tends to 37r/4, and the term cos 7 + sin 7 on 
left-hand side approaches zero. However, the ct" term on the left hand side grows faster, and 
therefore the entire left-hand side goes to 00. On the other hand, when a = 0, the left-hand 
side becomes 2v^/3. Therefore, varying ot the left-hand side obtains all the values in the 
interval (2-\/2/3, 00). This means that for any fixed e there exist infinitely many n values 
for which equation (12.151) has a solution a„. Consequently there are infinitely many real 
eigenvalues of order 0(5"^). 

The JBKW expansion above captured the real eigenvalues. They are all of 0[e~'^\ However, 
this expansion fails when a = and — 1 < /3 < 1 since in this case there is a turning point. 
We shall now construct a solution for such a case. If we order the eigenvalues by their real 
part, the eigenvalues we shall now construct come before those derived above. Typically a 
JBKW expansion captures the large eigenvalues and the associated oscillatory behavior of the 
eigenfunction. The lower eigenfunctions tend to oscillate less. Since we are now looking into 
the possible construction of a complex eigenvalue, we need to recall that they come in pairs 
of conjugate numbers. Geometrically it means that we anticipate two eigenfunctions, related 
by PT-symmetry, and therefore we anticipate in one case a turning point in [—1,0) and in 
another case a turning point in (0, 1], i.e. the turning point is at x = ±/5. While we seek now 
an eigenvalue that to leading order is purely imaginary, we shall be able to obtain also a lower 
order correction for it that will have real and imaginary parts. Therefore we assume that the 
eigenvalue is of the form 

(2.16) A~ie~2/?o + e"''(ai+«A), 
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Figure 2. The (numerically computed) real part of the leading eigenvalue 
(solid line) compared to the the asymptotic expansion (12.221 (dashed line). The 
horizontal axis in this figure is the current /, and the vertical axis is the real 
component of the eigenvalue. 

where the exponent v is still to be determined. 

Consider, then, without loss of generality, the case where the turning point is in the interior 
of the interval (—1,0), i.e. — 1 < /3o < 0. To leading order the potential is balanced exactly 
at the point x = —Po- We therefore expect the eigenfunction to be supported in a small 
neighborhood of this point, and to decay away from it. Hence, we construct an internal layer 
around this point by defining an inner variable y through 



largereal 



I73] (2.17) 



X = -/3o + e^y 



for some 7 > 0. Substituting the transformation (12.171) into equation (11.91) we see that in 
order to balance the different terms in the equation we must set 7 = 2/3 and u = 4/3. We 
thus obtain on the y scale the balanced equation 



"i75] (2.18) 



Uyy + iyU + (ttl + iPl) U = 0. 



While equation (I2.18P describes the internal layer form of the eigenfunction u, the outer 
solution is of course u = 0. Therefore, matching the inner and outer solutions implies that we 
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should consider equation fl2.18p over the entire real line with the conditions 

Iff] (2.19) u{y = ±00) =0. 

The eigenvalue problem (I2.18I) - (I2.19I) can be solved explicitly in terms of Bessel (or Hankel) 
functions. This was done by Ivlev and Kopnin [8] who concluded that this problem has no 
solution with finite L2 norm. This means that our assumption that /3o < 1 is not consistent. 
In other words, the concentration cannot occur in the interior of [—1, 1]. We therefore consider 
now the last remaining case in which jSo = 1. The scaling (12.171) of the internal variable y 
is the same, except that now it is more appropriate to call it the 'boundary layer' variable. 
Thus, the boundary layer equation (I2.18P is considered over the half line y G [0,cxd). The 
Dirichlet condition and the matching to the outer solution together imply the condition 

"i79] (2.20) u{0) = u{oo) = 0. 

The half-line eigenvalue problem (12.181) . (12.201) was also studied in [8]. In this case the eigen- 
function has finite L2 norm. The authors computed the leading eigenvalue to be approximately 

Isi] (2.21) ai + = 1.17-2.022. 

Returning to the original notation for the current, we obtain the eigenvalue 

l83] (2.22) A ~ 1.17/2/^ + - 2.02/2/3) for / > 1. 

In Figure [2] we depict the asymptotic expansions for the real part of the leading eigenvalue 
(dashed line) and the numerically computed real part (solid lines). Similarly, we depict in 
Figure [3] the asymptotic imaginary part (dashed line) and the actual imaginary part (solid 
line). In both cases the curves are very close to each other (the error is roughly 0.02). 

The eigenvalues are arranged by their real part in an ascending order. It is found numerically 
that as / increases towards the critical value of Jc ~ 2.27 (for L = 1 and Neumann boundary 
conditions) and at the critical value Ic ~ 12.31 (for L = 1 and Dirichlet boundary conditions), 
the first and second eigenvalue collide. These results are consistent with the bounds of ref. 
[12] . In the Neumann case their estimate for / below which the entire spectrum is real is 
/ < 71^/8, while the corresponding Dirichlet estimate is / < 3tt'^/8. 
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Figure 3. The (numerically computed) imaginary part of the leading eigen- 
value (solid line) compared to the asymptotic expansion (12.221 (dashed line). 
The horizontal axis in this figure is the current /, and the vertical axis is the 
imaginary component of the eigenvalue. largeimag 



3. Eigenvalue collision in PT-symmetric problems 

In this section we look in some detail at the collision process of two eigenvalues. Our goal is 
to derive an asymptotic expansion for the eigenvalues and eigenf unctions near collision, using 
/ — Jc as our small parameter. The analysis in this section will again be formal. 

For this purpose let / be near a critical value Jc the first two eigenvalues (ordered by their 
real parts) Ai and A2 collide; that is, the eigenvalues are real for / just below Jc but coincide 
at / = Jc. Let the associated eigenfunctions be ui and U2, respectively. We should note that 
the analysis we will present below is valid near the collision of any two eigenvalues but we 
focus on the first collision since this is most relevant to the stable bifurcation picture to be 
presented subsequently. 

It follows from equation (11. 9p that they satisfy 

s9 (3.1) J Ui{x)u2{x) dx = 0. 
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We term this property as PT — orthogonality. It is different, of course, from usual ortliog- 
onality. As long as the eigenvalues are real, their associated eigenfunctions are (up to a 
normalization) PT-symmetric. Therefore PT-orthogonality is the same as orthogonality in 
the Krein inner product: [f,g] = J f{x)g*{—x) dx = [12]. 

Now, as / approaches Jc, we assume that Aj and Xu approach a common value A(°). At the 
same time, Ui and U2 also approach a common function that we call u^^\ Notice, that this 
statement is made after some proper normalization, since the problem is linear. We really 
mean that the ratio Ui/u2 approaches a complex constant. To see why, assume to the contrary 
that there are two independent eigenfunctions Ui and U2 associated with the real eigenvalue 
A'^°^ Equation (11.91) is of second order, and therefore its solution space is spanned by two 
independent functions. However, ui and U2 cannot form such a basis, since they both satisfy 
homogeneous Dirichlet boundary condition, while clearly there are solutions of ( II. 9p that do 
not satisfy such conditions. 

Having established that the collision eigenvalue A*^°^ has an algebraic multiplicity 2, but a 
geometric multiplicity 1 (and thus we can say that the operator M is Jordan at Jc), we proceed 
to study what happens when I increases past I^- 

We first write the eigenvalue problem for the operator M near in the form 

sll (3.2) Uxx + i{Ic + £a)xu + Xu = 0, 

where the sign of the parameter a determines if we move up or down from Ic- Because of the 
singular nature of M at 1^ it turns out that the perturbation scheme is not analytic. Rather, 
we need to expand the eigenfunction u and the eigenvalue A in powers of e^^^: 

[il^ (3.3) m = mW + £1/Vi)+£m(2) + ....^ A = AW + £1/2^«^^^(2)^____ 

At the first order we find of course 
[il^ (3.4) := + ^x/cm(°) + A(°)m(°) = 0, 
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where we used this opportunity to introduce the operator notation L. It is important to note 
that the PT-orthogonahty noted above imphes 

(3.5) y rfx = 0. 



At the 0{e^l'^) level we get 
^ (3.6) /:«(i)=w« + zWi) + A(%« = -A«n(°). 



s21 



19 



21 



s23 



Multiplying the last equation by v!^^^ ^ integrating over the interval [—1, 1], and using the PT- 
orthogonality (13. 5p we see that equation (13. 6p is solvable. However, unlike the case of regular 
eigenvalue perturbation schemes, we gain no information on A'-^-' at this level. We therefore 
need to proceed to the 0[e) level: 

(3.7) ^ ^(2) ^ ^^^^^(2) ^ ^(0)^(2) ^ _^(1)^(1) _ ^(2)^(0) 

To get a solvability condition we multiply both sides by ■u^''^ integrate over the interval and 
use (13. 5p to find 

(3.8) "^'^(^)"^°^(^) = ^"/ a;M(°)(x)2rfx 
It is convenient at this point to introduce some notation. First, we set 

(3.9) = ReM(°)(x) + iImM(°)(x). 

Using the PT symmetry of n*^'^-*, we choose a normalization in which ReM*^°^ is even, while 

Ijn-yCo) is odd. Then we define the real parameter a\ through 

1 /.I 

1 



(3.10) ai = -ij xu^^\xf dx = 2 j xReu^^\x)\mu^^\x) dx 

Next, let K{x) be the solution of the nonhomogeneous ODE 
^25] (3.11) CK = K^^ + ihxK + \^^^K = u^^\ K{±l)=Q. 



The identity (13. 5p ensures that equation (13. lip is solvable. Using this canonical function K, 



we express u^^'^ as 



^rf] (3.12) mW(x) = -\^^^K{x). 
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Finally, we define 

"i29] (3.13) b = j K{x)u^^\x) dx. 

Notice that K is also PT-symmetric, i.e. K{x) = K{—x). In particular Rei^' is even while 
Im K is odd. A numerical integration of K, and a numerical evaluation of the functionals in 
(KW) and (137[3|1 gives 

l30] (3.14) ai^O.29, 6^0.12. 

We use the notation above to derive from equation (13.61) the relation 

Jsi] (3.15) = -aai/b ^ -2.42a. 

When a < 0, J is just below Ic and there are two real solutions. The negative one corresponds 
to the first eigenvalue, and the positive one corresponds to the second eigenvalue. On the 
other hand, when / increases past Ic, i.e. when a > 0, there is a complex pair of conjugate 
solutions. This implies that the critical eigenvalue A'-"^ splits into a complex conjugate pair 
with the 0(6^^"^) correction A^^-* being purely imaginary. We remark that the same analysis 
applies to any collision of real eigenvalues. In Figure H] we compare the asymptotic expansion 
(13. 3p of the imaginary part of the first eigenvalue (dashed line) with the numerically computed 
value (solid line). 

To find the next term A^^-* in the eigenvalue expansion, we proceed further to the 0(£:^/^) 
level: 

(3.16) £^(3) = -A(2)m(i) - A(i)^i(2) _ ;^{3)^{o) _ ^^^^(1)^ 

Multiply equation (I3.16P by u^^'^ and integrate by parts over [—1,1] to get 
l63] (3.17) Aa / u^^'^u^^Ux + X^^^ [ u'^^^ u^^Ux + za [ xu^^^u^^Ux = 0. 



J~i J-i J-i 

Substituting the relation (13.121) into equation (13.171) and dividing by A*^^^ gives 

165] (3.18) -bX^^^ + J u'-^^u'-^Ux-ta J xK{x)u'-°Ux = 0. 
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0.1 ' ' • • • ' • • ' ' 

12.32 12.33 12.34 12.35 12.36 12.37 12.38 12.39 12.4 12.41 



Figure 4. The (numerically computed) imaginary part of the leading eigen- 
value (solid line) compared to the leading order term in the asymptotic expan- 
sion (13.31 (dashed line). The horizontal axis in this figure is the current /, and 
the vertical axis is the imaginary component of the eigenvalue. 

It is useful at this point to introduce additional canonical functions and functionals in the 
spirit of K, ai and b defined above. Thus we define two canonical functions ( and w through: 



M3] (3.19) CC = K{x)-ei. c(±i) = o. 

Here 6i is a constant chosen such that equation (13.191) is solvable. Namely 

1 .1 
K{x)u^^\x)dx/ / u^^\x)dx. 



h43b 



-1 



(3.20) Oi 
The next canonical function w{x) is defined by 
his] (3.21) Cw = -ixu^^^ - 02, w{±l) = 0, 



h45b 



with 
(3.22) 



62 = -I xu^^\xfdx/ u^^\x)dx. 



-1 



Using (, w, and K we can write 
^ (3.23) u^^\x) = (X^^^yCix) + aw{x) + X^^^K{x). 



image lam 
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We further define the functionals 
h45c| (3.24) ^ j ^"''"^ ^'^^ j d.3 = ij xK{x)u^^\x) dx, 

A numerical computation gives 

h45d] (3.25) di ^ -0.014, d2 ^ -0.02, d^ ^ -0.02. 

Using these functionals and formula (13.231) for u^"^^ in equation (13.181) . and using (I3.15P to 
eliminate the contribution of the coefficients 61,62, gives the following expression for A2: 

Ist] (3.26) A^^) = (ads - 0^2 + {X^^^)^di) /2b. 

We can conclude now an interesting fact. The eigenvalue A*-'^-' splits into two eigenvalues as 
the current is varied away from Ic- For I < Ic^we obtain a real pair, while for I > Ic^e obtain 
a complex pair. This splitting manifests itself in the two values for A*^^-* obtained from equation 
(]3.15p . We can also see how the single eigenfunction u^^^ splits into two eigenf unctions through 
equation (13.61) . On the other hand, A''^^ appears in equation (I3.26P only through its square. 
Therefore equation (13.261) implies that A^^-* is unique and real. In particular, if we draw the 
real part of the colliding eigenvalues as a function of / near the collision, we obtain that the 
function is not analytic at /~. In fact dXo/dl blows up as we approach Ic from below, due to 
the 0{e^^'^) contribution to the expansion (13.31) coming from (13.151) . Yet due to the fact that 
A*^^^ is purely imaginary for / just above Ic, we see that the graph of Re Ai(/) (= Re A2(/)) is 
differentiable from the right at / = Ic- This analytical conclusion is verified in the numerical 
solution. 

In Figure [5] we compare the asymptotic expansion (13.31) of the real component of the first 
eigenvalue (dashed line) with the numerically computed value (solid line). The two lines are 
almost indistinguishable, and the error is 0(0.001). 



4. The bifurcation from (N) to (S) 

In the next two sections we study the shape and the stability of the bifurcation branch from 
the normal (N) state to either the stationary (S) or periodic (P) state using formal asymptotic 
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Figure 5. The (numerically computed) real part of the leading eigenvalue 
(solid line) compared to the leading order term in the asymptotic expansion 
(13.31 (dashed line). The horizontal axis in this figure is the current /, and the 
vertical axis is the real component of the eigenvalue. reallam 

expansions and multiple time-scales. Later, in Section 6, we will present a rigorous justification 
of these calculations by appealing to center manifold theory. 

For this purpose, it is convenient to rewrite the system f ll.4l) -( |T31) as a single nonlocal 
complex equation by first solving for the electric potential in (II. 5p . thereby obtaining 



i r 

solvephi (4.1) = -Ix + - (V'V'x ~ V'*^x) dx' . 
2 Jo 

Then we substitute this into (11. 4p to obtain 

(4.2) i/jt = iIj,, + ixi^j + rv- + U[^i 

where we have introduced notation for the cubic nonlinearity 

(4.3) M[i^] := - ^ + i ^ r (^^: - ^>,) dx'. 



ref orma 



nonnot 



As always, this equation is augmented with Dirichlet boundary conditions at x = ±1 and 
initial conditions. 



Lone 



ref ormb 



h279 



h281 



h285 
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We work in this section with a fixed current I in the regime I < 1^ where the leading 
eigenvalue Ai of the operator M is real (as is the entire spectrum) and we expect a bifurcation 
to a stationary state. The transition takes place exactly when F crosses the value Ai. This 
defines the curve V = ri(/). 

We denote the leading eigenfunction by Ui, and normalize it as usual by Ui{0) = 1. Recall 
that Ml satisfies the equation 



(4.4) 



Liui 



+ ixlui + XiUi = 0. 



To find the solution to (14.21) just above the curve ri(/), we set F = Ai + e, so that (14.21) 
takes the form 



(4.5) 



Anticipating the contribution of the nonlinear terms in the forthcoming expansion, we seek 
a solution that is proportional to leading order to ui plus a small perturbation: 



(4.6) 



iIj{x, t) ~ £:^/^a(r)ui(a;) + e^^'^tpiix, r) + 



Since F is a small perturbation of Ai, we expect the time evolution to be slow, hence we have 
introduced the time-scale r = et. Our goal here is to compute the function a(r) and thus to 
obtain completely the leading order term in the expansion. 

Substituting the ansatz (14. 6 p into (14.51) . we see through (14.41) that the 0(e^/^) terms are 
balanced by the choice above for ip, with the function a{T) not yet determined. Proceeding 
the 0(6^^"^) level, we obtain 



(4.7) 



-ar + a)ui+ N[oiUi\. 



To obtain a solvability condition for ipi, we multiply equation (14.71) by Ui{x) and integrate 
over [—1, 1]. We obtain the following equation for a(r): 

/^-^ N'[aui\ui dx 



(4.8 



a + 



a + a. 



h287| (4.9) 



j^^ u\ dx 

After a lengthy calculation, the coefficient xii is found to be given by 

1 



-ciiii - 7ii //5 
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where 
h289] (4.10) 

/I /*! pi px 

u\dx, 7ii = / \ui\^u\, Clin = / ulOu dx, 9ii = / Ui(mi)* - ^^(mi)^ dx'. 
1 J-i J-i Jo 

The stabihty of the bifurcation branch depends on the sign of Xii- Examination of (14.81) 
reveals that there is stable branch of equilibria if Xu < 0) with |a| = ^^j^ . We note that in 
light of the rotational invariance of the whole problem (14. 2p . there is in fact an entire circle 
of equilibria with a given by ^ — e'^°, 6*0 G [0, 27r). 

On the other hand, if xii > then an unstable branch of equilibria exists for e small and 
negative wherein (14. 8 p is replaced by 

Ot- = —a + 

and the equilibrium value of a is given by l/^/xn- It turns out that both signs can occur, 
depending on the current /. For instance, when / = 7, we get 

(3(7) = 0.785, 7ii(7) = 0.652, ciiii(7) = 0.375, Xii(7) = -0.592, 

while 

= 0.403, 7ii(ll) = 0.449, Ciiii(ll) = 0.924, Xii(ll) = 0.03. 

A careful computation shows that xii vanishes at / ~ 10.93. We denote this critical value 
by Jfe and conclude that the bifurcation from (N) to (S) is stable (type II) for I < 1^ and 
unstable for Ik < I < Ic- Finally we point out that the expansion in this section breaks down 
for / near J^, and one needs to proceed to higher order terms there. 

5. The bifurcation from (N) to (P) 

In this section we compute asymptotic approximations for the solution in the oscillatory (P) 
state. In this state the solution ilj{x,t) is time-periodic. We assume in this section that the 
current / is fixed in the regime I > Ic- The transition to the (P) state takes the form of a Hopf 
bifurcation; namely, the real part of the spectrum is zero, and the bottom of the spectrum 
consists of a conjugate pair of purely imaginary eigenvalues. We point out that the spectrum 
of the operator M (cf. (11.90 ) does have a nonzero real part, but this real part is exactly 



NP 
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balanced at the transition curve by our choice in this section of F = ri{I)+e = KeXi + e with 
£ > in (14. 2p . Hence, if we extend our definition fl4.4p of the hnear operator Li to include 
the case where Ai is complex via 

(5.1) Liu := Uxx + ixlu + (Re Xi)u, 

then it is Li that possesses a pair of purely imaginary eigenvalues with corresponding eigen- 
functions Ui and U2 satisfying 

imeig (5.2) LiUi = —i(lm \i)ui and L1M2 = ^(Im Ai)m2. 

We choose to normalize the eigenfunctions so that Uj{0) = 1, j = 1,2, and we assume Ai is 
defined so that Im Ai > 0. With the above choice of F and definition of Li, we again find that 
dllD-dLSl) takes the form M . 

At leading order, we expect the solution to be comprised of a linear combination of solutions 
to the equation ipt = Liip. This leads us to seek a periodic solution to (14.51) of the form 

h79] (5.3) ^{x,t) e'/^aiiT)e-~'^'^^''ui{x) + a2iT)e'^'^^''u2ix)) + e^/^ijiix,t) + ... 

Here ipi is assumed to be a periodic function of t with period = 27r/ImAi + 0{e) and 
ai,a2 are coefficients that we expect to evolve slowly in time. We thus have set ai = Q;i(r), 
where, just as before, r = et. Our goal here is to compute the functions ai(r) and thus to 
obtain completely the leading order term in the expansion. 

Substituting the ansatz (15. 3p into equations (14.50 shows that the 0{e^^'^) terms are balanced 
by the choice above for ijj, with the coefficients aj(r) not yet determined. Proceeding to the 
0(5^/^) level, we obtain 

(5.4) {^,)t-L^iJ,= e-'^'"^iV(ai-«ir) + e'''"^^*M2(a2-a2r) 

hsTI (5.5) + Ar[aie-^^'^^i*Mi + a2e^^™^^*M2], 



cf. (I43D. 

To obtain a first solvability condition for ipi, we multiply equation ( 15. 5p by e*^™'*'^*'Ui(x) and 
integrate over [—1, 1] x [0, 27r/Im Ai]. A second solvability condition is obtained by integrating 
similarly against the function e~*^™'^^*M2(a;). We note that through (15.20 and the assumed 
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periodicity of ipi, the first integration against the left-hand side of fl5.5p yields, after an 
integration by parts: 

»1 /•27r/ImAi 



/I />27r/ImAi rl /•27r/Im Ai 

/ tPie'^'^^^^uidtdx- I j iJie^^'^^^^LiUidtdx = 0{e) 

1 Jo J~i Jo 



Similarly, the left-hand side in the second integration vanishes to leading order. 

After a lengthy but straight-forward calculation, these two integrations then give rise to a 
pair of equations for the coefficients ^^(t), namely 

h85a| (5.6) {ai)r = ai + {xii\ai\^ + Xi2\a2\^) ai, 

h85b| (5.7) (a2)T = "2 + (Xii|«2|^ + Xi2|ain «2- 

Here the coefficient xii is again given by (14.91) while Xi2 is defined by through: 

Ii87| (5.8) Xi2 = ( Jcii22 + Jci2i2 - 2721 ) //?, 



2 2 



where 
h89] (5.9) 



1 /•! rl 



(3 := u\dx, 'jij = / \ui\^u]dx, Cijki = / UiUjOki dx, 9ki = / {uk{ui)l - u*i{uk)x) dx' . 



1 J-l Jo 



We note that the notation above is consistent with ( ]4.10p . 

To analyze the evolution of the a-s, we note that from fl5.6l) - fl5.7p it is easy to derive the 
system 



mh85a| (5.10) (|ai|)^ = |ai|(l + (Rexii|«ir + Rexi2|a2r), 

mh85b | (5.11) (|a2|)^ = |a2|(l + (Rexii|a2r + Rexi2|ain, 

governing the evolution of the moduli of the a^. Now if the initial conditions for (14.51) are 
taken such that ai(0) 7^ and ^2(0) 7^ 0, then it follows from fl5.10l) - fl5.11l) that ai and 0:2 are 
non-zero for all future times. In this case, we introduce the Ricatti transform r := |ai|/|a2|- 
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From flCT?]) - fl3TT|) we obtain that 



r 







/ 


"2 


/ 






^2 


2 



req 



(5.12) 





|a2|(Rexii - 


Re 


Xi2)(|ai 


2 _ 


«2 






a2 


2 



= (Rexii - Rexi2)\a2\Mr^ " 1) = "(Rex) |a2rr(r2 - 1), 
where denotes ^ and we have introduced the complex constant 



chihatdefnl (5.13) X '■= Xi2 - Xu- 



Provided that Re^ > 0, we learn from (15.121) that r ^ 1 at an exponential rate as r — > oo. 
One can indeed check numerically that for / > Jc, the inequality Re x > holds. See Figure 

El 

Returning to the system f l5.10p -( l5.1ip with this information, we can determine the asymp- 
totic value of the modulus of both ai and 0:2 to be 



|ai(r) I ~ |a2(^) I ~ Wtf; ^Ti \ = \ TT—:^ for r >> 1 

^' (Rexii + Rexi2) V Rex 

where we have introduced another complex constant 



itildedefn] (5.14) ^ ■= -{Xn + Xu)- 



Numerical calculation reveals that Re x > for I > 1^ as well. Again, see Figure El 

Substitution into fl5.6l) - fl5.7p then yields a linear dependence on r of the phase of both ai 
and a2 for r large with 



asymal I (5.15) aAr) = alir) = J e'"^^ for r >> 1, 

Rex 



where w := Of course, in light of the rotational invariance of all of the above equations, 

the asymptotic state (I5.15P holds only up to a constant rotation. As an example, we computed 
the functionals xii ^"^^ X12 for / = 20. In this case the leading eigenvalues are Ai^2 ~ 
8.64 ± 5.25 i, the amplitude is ~ 0.92 and uj ~ —1.81. Unlike the case I < Ic, the transition 
from (N) to (P) is always stable. 
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Figure 6. The dashed curve represents the graph of the parameter Rex as a 

function of apphed current / and the sohd curve represents the graph of Rex- chi_l 

The case where either a;i(0) = or 0:2(0) = is treated separately. In view of fl5.6l) - fl5.7p . 
if, for example, ai(0) = 0, then necessarily Q;i(r) = 0. Hence, 02 evolves by 

{a2)r = (1 + Xll|tt2p) "2- 

From this and (15. lip it easily follows that 

^ Imxil 

a2unst | (5.16) 0^2 (t) ~ - eR<=xii ^ for r >> 1 

V-Rexii 

provided that Rexii < 0. Since xu = ~^{x + x); aiid both x X have been shown 
numerically to be positive, we do indeed have this condition on xu met. Similarly, when 
^2(0) = 0, one finds that a2{r) = and 

]L _ ■ r™xii 

alunst (5.17) '^ij'^) ^ - p'^Roxh forr>>l 

V -Re Xii 

provided again that Rexii < 0. 

To summarize, in the generic case where neither ai(0) nor 0:2(0) vanish, the solution in the 
(P) state is to leading order 

h95] (5.18) ^{x, t) ~ e^/^A (e-'(^'"^i+"^)*Mi(x) + e''-''^^'+^'^'u2{x)) . 
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Therefore the solution is time-periodic with period 



95l] (5.19) ~ 27r/ (Im Ai + ue) + o{e) 



In a subsequent article, we will examine in some detail how the asymptotic solution (15.181) 
extends deeper into the nonlinear regime of the (P) state. 

6. Rigorous bifurcation theory 

In this section we will establish a rigorous justification for the formally derived expansions 
and bifurcations of the previous two sections. As was done earlier, we will set the normal 
conductivity a equal to one and set the interval [—L,L] to be [—1,1], thus focusing on the 
interplay between the externally forced current I and the temperature-dependent parameter 
r. As was done earlier, we will treat the case of Dirichlet boundary conditions, though 
similar conclusions can be rigorously established for the case of Neumann boundary conditions. 
Though in our previously derived formal asymptotics, we only pursued the case of bifurcation 
off of the principal eigenvalue, at the end of this section we also treat bifurcation off of any 
eigenvalue. 

6.1. Spectral properties of linear operator. We begin by recalling some notation and 
collecting some facts about the linear operator M given by (II. 9p . 



Mspec Lemma 6.1. The spectrum of M consists only of point spectrum, denoted by {— Aj} with 
corresponding eigenf unctions {uj}. If{Xj,Uj) is an eigenpair satisfying 

Muj = —XjUj, ^^(±1) = 0, 

then 



Relm (6.1) ReAj>0, and |ImAj|</. 



Thus, in particular we may order the eigenvalues Ai, A2, . . . according to the size of their real 
part, with < Re Ai < Re A2 < . . . . The PT-symmetry of the operator is reflected in the fact 
that if {Xj,Uj) is an eigenpair then so is (A*,Mj) where uUx) := u*{—x). Finally, for each 
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positive integer i, there exists a positive value of the current I , which we denote by Ii, with 
li < Jf+i such that 



split I (6.2) X2£-i, >^2i e K for I <h while \2i-i = >^li ^ M for I > h 



Remark 6.2. We note that in earlier parts of this paper, the critical value Ji was denoted by 
Jc- In the next subsection, we revert to this notation to keep consistency with earlier sections. 

Proof. The fact that the spectrum consists entirely of eigenvalues follows from standard spec- 
tral theory. To verify (16. ip . multiply the equation Mu = —Xu by u* and integrate to obtain 

jli\uxf dx I J^-^^x\uf dx 
ReA = — I ^ and ImA = j ^ . 

\u\ dx 1^1 dx 

The fact that if (A, u) is an eigenpair then so is (A*,m^) follows by inspection. The assertion 
that the spectrum is real for small / follows from Theorem 4.1 of |12| while the existence of a 
critical values {le} beyond which pairs of real eigenvalues collide to form conjugate pairs is a 
result of [13115]. □ 

Next, for any fixed positive integer i, we introduce the operator 



Lell| (6.3) LeU := Mu + (Re Aa^-i) u. 



Then Li has spectrum shifted from that of M by ReA2£-i so that L^Uj = —fijUj where 



Hj = \j — Re A2£-i. Based on the behavior of the spectrum of M described in Lemma [6. ![ we 



Iless 



eig 



gap 
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have the following scenario for L^. 

(6.4) For < / < /£ : 

/i2£-i = 0, Re/ij < for 1 < J < 2£ - 1, fij > for j > 2i - 1. 



Izj (6.5) For I = h: 

/^2^-i = /i2£ = 0, Re yUj < for 1 < j < 2£ - 1, > for j > 2£. 



Imore (6.6) For I > I^: 

Re/X2^_i = Re^.2t = 0, lmfi2i-i = -Im/i2^ ^ 0, 
Re < for 1 < J < 2£ - 1, Refij > for j > 2i. 

We also note that by PT-symmetry, if we normalize all eigenfunctions so that ^^(0) = 1, 
then we must have 



(6.7) U2e^i = 'u^2i-\ when I < le while U2e-i = ul^ when I > Ii. 
For later use, we also introduce the spectral gap (e given by 

(6.8) Q := inf {|Re fij\ : fij is an eigenvalue of Le with non-zero real part}. 

j 

In what follows it will be convenient to choose a basis for the eigenspace of Li that is 
PT-symmetric. This has already been taken care of when I < Ii. However, for / > Ji we 
introduce the basis Vi and V2 given by 



vdef n I (6.9) fi := Ui + U2 and V2 ■= i{ui — U2)- 



The PT-symmetry of this basis follows from (16. 7p . Note also that this basis satisfies the 
relations 



veval (6.10) LiVi = —Im X1V2, Lif2=ImAifi 



reform 



nonnot 1 



cubiclem 



cubic 
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6.2. Bifurcation from first eigenvalue. We now develop the rigorous bifurcation theory 
associated with the stationary and periodic solution branches formally derived in earlier sec- 
tions. To this end, we wish to reformulate the full nonlinear system f ll.3p -( fL6ll in such a way 
as to make it amenable to standard center manifold and bifurcation theory. Accordingly, we 
first solve for the electric potential in (11.51) . thereby obtaining 



solvephil (6.11 



(V^C - dx'. 



In this section, we will focus on the stable bifurcations that occurs off of the first eigenvalue 
of the linear operator and so we henceforth fix the positive integer i from the previous section 
to equal 1 and pick T to be of the form F = Re Ai + e. Substituting (14. II) into (11.41) . we obtain 
a single, nonlocal complex equation so that (ll.4l) - (ll.5p can be rewritten as 



(6.12) 



tljt = Ll^P+^f{^^J,e) 



where we have introduced Af{y, e) := Af[y] + ey with 



(6.13) 



^^[y] ■=-\y\ J ~ 



We recall that as before, the system is augmented with Dirichlet boundary conditions ip{±l,t) = 
along with the normalization 0(0, t) = for t > 0, and initial conditions, say iplx, 0) = ipo{x). 
We take e to be small and positive (unless otherwise specified). 

We remark that since Liy"^ = {Liy^ and A/'[?/^] = {Af[y])\ it follows easily that the flow 
(I6.12P preserves PT-symmetry in the sense that if t/' is a solution to (16.120 . then so is tp'' . 
Hence, by uniqueness, we note that if the initial data ipo is PT-symmetric, i.e. if ipo = ip^, 
then so is the resulting solution ip. 

For later use, we also record the estimate: 

Lemma 6.3. There exists a positive constant Cq such that 

(6.14) my]\\m<Co\\y\\%. for ally EH^{i-l,l)y,C). 

Proof. For y G Hq{{—1,1)) recall that both ||?/||^oo and \\y\\fji are controlled by ||?/a-||^2 . We 
begin by estimating the norm of the local part of A/". We find 



\y\ y\ 



L2 



< WyWle < \\y\\la. \\y\\L2 < C \\y\\\i . 



ezero 
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Similarly, 
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Turning now to the nonlocal part of M we find 



y / {yy*x-y*yx)) dx 



L2 



i\y\ \y\ \y^\ dx'Y dx 



1/2 



<C{j'_^ \y\ \yA dx) \\y\\^, < C Ml. \\y4^, < C \\y\\%, . 



Finally, we check that 



y / iyy*x - y*yx)) dx' 



< 



L2 



\y\ \yx\ dxfdx^ + ^(/i '^^ 
< c IblLa \\yx\\l2 + C \\y\\loo hxh^ < C \\y\\%i 



1/2 



□ 



Armed with a full understanding of the linearized operator and control on the nonlinear 
operator provided by Lemma 16.31 we can proceed to construct a center manifold for the flow 
fl6.12p . To this end, we will denote by Sc the center subspace associated with Li, that is, Sc is 
the eigenspace associated with any eigenvalues of Li having zero real part. For I < Ic {= h), 
we have Sc = spanjui} while for / > J^, Sc = span {t>i, t>2} (= span M2}). 

As e will play the role of a bifurcation parameter, we then augment (16.121) with the equation 



(6.15) 



et = 0. 



The theorem below provides for the existence of a finite dimensional invariant manifold as- 
sociated with the flow (I6.12p for each fixed small e describing all orbits of sufficiently small 
norm. This, in effect, allows us to rigorize the formal bifurcation calculations of the previous 
two chapters by reducing the analysis of the nonlocal P.D.E. to a study of a local system of 
O.D.E.'s with accompanying rigorous error estimates. 



cm I Theorem 6.4. For each value I > and positive integer k, there is a local center manifold 
Ai C H^{{—1, 1); C) X M of (I6.12p . (16.150 tangent to the center subspace. The center manifold 
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Ai is expressible as a graph over the center subspace in the sense that there exists a map 
$ : 5c X M ^ ^^((-1, 1); C) such that 



cm; 



|] (6.16) = [J {Me^{e}) where Me ■= {^{u,e) : u e S^, \\u\\^^ < 5q, \e\ < Sq] 

|e|<£0 

for some sufficiently small positive constants 5q and Eq depending in particular on k. The 
center manifold is invariant under complex rotation, i.e. {ip,e) G A4 =^ {e^^ip.e) G for 
all 6 eM. and in fact 



phirot (6.17) e'"^{u,e) = ^{e''u,e) forallOeR. 



The center manifold is also PT- symmetric, i.e. {ip,e) G Ai =^ {4'^^) ^ Ml. If u = 
then ^{u,e) = <l>"^ («,£:). 

The discrepancy between the center manifold and the center subspace can be expressed 
through the estimate 



very tang I (6.18) ||$(m, £:) - < Ci + |£:| ||n||^i) 



which holds for any pair {u,e) such that u E Sc with \\u\\fji < 6o and \6\ < Eq, where Ci is a 
positive constant independent of u and e. 

The center manifold is locally invariant for the flow (16.121) in the sense that if \£\ < sq and 
the initial data ipQ lies on Me, then so does the solution ip'^ to (16.121) so long as \\'ip'^{-,t)\\jji 
stays sufficiently small. Hence, for such initial data, one can describe the resulting solution 
ip^{t) = ip'^{-.,t) through either one or two maps /?? : [0, oo) — > C via ip^{t) = $(/3^(t)-Ui, e) 
when I < Ic or tp^{t) = $(/5f (t)fi, /3|(t)f2, £:) when I > 1^). Finally, M contains all nearby 
bounded solutions of (16.121) in , and in particular, it contains any nearby steady-state or 
time-periodic solutions. 

Proof. We follow a standard center manifold construction, along the lines for example, of [3]. 
To outline this approach, we first note that in light of conditions (16. ip . the spectrum of the 
operator —Li lies within the set 

{AgC: |arg(A + a)|<^}, 
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for some positive number a = a{I). Hence, —Li is sectorial and we may assert the existence 
of an analytic semi-group {e^^^}t>o, cf. [7j, Theorem 1.3.4 or JT9\, section 2.2.3. 

We will denote the (L^) projection operators from H^{{—1, 1); C) onto the center and stable 
subspaces of Li by and lis respectively. Here by stable subspace we mean the span of all 
eigenvectors of Li whose corresonding eigenvalues have negative real aprts. We note that since 
the real part of all eigenvalues of Li are non-positive, Li has no unstable subspace. 

A local center manifold is constructed by first constructing a global center manifold for 
a problem with a truncated nonlinearity through the introduction of a cut-off function p G 
C°°{[0, oo); [0, 1]) satisfying p{s) = 1 for < s < 1 and p{s) = for s > 2. For any 6 > we 
then let p^{s) := p{s/6). We use this cut-off to truncate the nonlinearity Af{y,e) by defining 
Af'{y,e) ■.= p'{\\y\\^,)Af{y,e). 

The graph map $ : x M ^ H^[[—l, 1); C) is then defined by the following procedure: 
For u & Sc and e fixed, we use a "variation of constants" approach, rephrasing the P.D.E. 
fl6.12p (with the original nonlinearity replaced by Af^{y,e)) as an integral equation: 

y{x,t) = T{u,e,y) := e^^^u+ [ e^'^^^'^^U^Af^ {y{x,T),e) dr 

Jo 

+ 1 e^i(*-")H,A/'*(y(x,r),e) c/r. 

J —oo 

(6.19) 

Then one argues that there exists a unique fixed point = yu,^{x,t) to fl6.19p in the space 
of functions that grow sufficiently slowly at t = ±oo given by 

Yr, := {y e C {-oo,oo)-H\{-l,l))) : < oo}. 

Here \\y\\^ ■= sup^gig e"'''*' ||z/(', ^)|Ihi((_i i)) and r] is any fixed positive number less than the 
spectral gap Ci, cf. ( 16. Sp . Once the existence of this fixed point is established, we define the 
map $ by 

We should remark that when t is negative, we interpret e^^*Hc in f l6.19p to mean fiow 
projected onto the finite dimensional center subspace; thus it reduces to a finite number of 
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ordinary differential equations. To see this and to carry out the apphcation of the contraction 
mapping principle to fl6.19p . one considers the inverse Laplace transform representations 

Mil (6.20) e*-^in, := [ e^\\I - Li)~^ dX, 

and 

bds] (6.21) e*^^n, := [ e^\XI - Li)'^ d\, 

where Tc is a bounded contour enclosing the eigenvalues with zero real part and Tg is a contour 
in the left half-plane enclosing the stable spectrum that tends asymptotically to infinity along 
the lines a(l — s) ± ais as s ^ ±oo. Since the resolvent is bounded along these contours, it 
follows from these representations that for every si G (0,Ci) one has 

"bd4] (6.22) ||e*^inc||^i^^i < C, for all t e M, 

bds] (6.23) ||e*^^n,||^,_^, < Ce~(^^~"i)* for all t > 0. 

Invoking these bounds, and by choosing the parameter 6 in the cut-off of the nonlinearity 
sufficiently small, the existence of a fixed point to (16.191) follows from (I6.14p by the contraction 
mapping principle, from which we find easily also Lipshitz regularity of $. The asserted C°° 
regularity of $ may be established by a careful iterative argument as described in [3], using 
(jco rgguiarity of the truncated equations (in general, the center manifold inherits one degree 
less regularity than the underlying equations); we omit discussion of this delicate point. The 
local center manifold for the untruncated problem is then realized through (16.161) by choosing 
6q and Eq sufficiently small. 

The rotational and PT invariance of A4 follow from the fact that the center subspace Sc 
enjoys these invariances and the fact that for any £ one has 

T{e''^°u,e,e''^°y) = e'''°T{u,€,y), as well as T{ul€,y^) = {T{u,e,y))^ . 

With regard to this last assertion, note in particular that N^{y\e) = [Af^{y,e)y . Also if 
u = m"!" and if yu,e = ^{u, e, yu,e) then necessarily 

vie = (r(^, e, yu,e)y = r(M, e, yl ) 
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and so by the uniqueness of the fixed point, necessarily yl^ = yu,e- Hence, in particular 
2/i,e('5 0) = Z/«,e(-5 0) and we have = Similarly, e*^o$(M,e) = $(e*^ou,£) for all 

Oo e M. 

Finally, we turn to the verification of fl6.18p . This comes from an examination of the 
iteration procedure leading to the fixed point as follows. Picking 6q sufficiently small, we may 
argue that, for instance, 

\\r{u,e,yi) - r{u,e,y2)\l^ < \ llz/i -1/2!!^ 

for all u G Sc, all sufficiently small e and all yi, y2 ^ Yri- Letting ?/„ denote the solution to the 
linear problem, i.e. Uu '■= e^^^u, it then easily follows that 

\\yu - yu,e\\^ = \\yu - T{u,e,yu,e)\\^ < \\yu - T{u,6,yu)\\^ + \\T{u,6,yu) - T{u,e,yu,e)\\^ 
and so 

hu - yuA\n ^ 2 \\yu - r(u,£:, . 

Then we calculate 

- u\\jj^ = \\yu,ei-,0) -?/„(-, 0)||j^i 

< sup -yui-,t)\\jji e~''l*l = \\yu,e - yuL 
teR 

< 2 \\yu - T{u,e,yu)\\^ < 2C (||n||^i + \e\ ||u||^i) , 
where in the last estimate we invoked (I6.14p . (16.191) and fl6.20l) - (l6.21l) . 

□ 

We also will need a version of the standard result on exponential attraction to an orbit on 
the center manifold in the absence of any unstable manifold. Again the proof we sketch is an 
adaptation of a more general but somewhat weaker result in [3] that is valid in the presence 
of an unstable manifold. 



exp I Theorem 6.5. For any positive integer k and r = r{k) > sufficiently small, there exists a 
map from B{0,r) C to Ais; equal to the identity when restricted to Aie, such that, 
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for all solutions tp^ of f l6.12p originating at time t = within B{0,r), 
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(6.24) ll^-(t) < Cie-n'dHi{r{0),Me 



so long as ip^ remains in B{0,r), where ijj^ G Aig denotes the trajectory along originating 
at time t = at Pi;{iIj'^ (0)) and f] > and Ci are uniform constants. Here dn^i-, ■) denotes 
the distance in H^. 

Proof. As noted in [TB] , this follows by the proof of the more general approximation property 
(v) of the Center Manifold Theorem stated in [3], restricted to the case that the underlying 
linearized operator (Li in this case) has no unstable manifold. 

To say a bit more about our adaptation of the approach presented in j3j, given a solution 
ip^ = ■ip'^[x,t) to (16.121) . one first extends ip'^ to a function ip^ defined for negative t- values via 



ijj^ for t > 
ip%Q) for t < 0. 



Thus, 4"^ represents a globally bounded solution to the equation = Liip^ + J\f {■ip'^ , e) + 
where 



phiL (6.25) (p' = (p'{x,t) :-- 



for t > 

-Li^^(O) -7V(^%£) fort<0. 



Then we fix any positive rj such that rj < (i (cf. (16.81) ) and seek a function z such that 
ip^ + z & Ai^. We will find such a 2; in the set 

Zr, ■.= {zeC ((00, oo);H\{-l, 1))) : |z|^ < 00} 

where |z|^ := sup^^^e'^^ \\z{-,t)\\j^if^i^__^-^^^-^ and then define the projection onto Ais via 
Psl'ip'^iO)) := tp'^lO) +z{0). Thus, the trajectory on the center manifold satisfying (16.241) will be 
ip^(t) := ip^{t) + z{t). The function z is produced as follows: plugging ij)^ + z into the integral 



newf ixed 



zeieq 



Flamb 
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equation (16.191) leads one to seek z as a fixed point of the mapping A : Z^i ^ defined by 

/oo 
gLi(t-r)n^ ^^S ^) + z{x, t),6)- t) , e)] dr + 



e^i(*-")n, (f)%x, T)dr+ / e^i(*-")n, [Af^ (^/'^(x, r) + z{x, T),e) - Af^ {ij'{x, r), e)] 



oo 



(6.26) 

Again the existence of a (unique such) fixed point z & follows readily from the contraction 
mapping principle since one can check that |A(2;i) — A(2;2)|r, < 6*1^1 — Z2\ri for some 9 G (0, 1). 
The fixed point z, and thus the map Pe that it determines, is Lipshitz in ip'^lO) by construction. 
With further effort, it may be shown to be for any k, by a procedure similar to that used 
to show smoothness of the center manifold [3] in the analogous fixed-point construction of 
Proposition 16.41 using regularity of both the center manifold and the truncated equations. 
(In general, inherits the regularity of the center manifold.) Note that the center manifold 
solution tp'^it) so constructed satisfies the truncated equations (16.191) and not (I6.12p . since 
the righthand side of (I6.26P involves the truncated nonlinearity Af^ in place of Af. However, 
this makes no difference since the equations agree on the ball 5(0, r) under consideration, for 
r > sufficiently small. 
Hence, we have 

\zU < k - A(0)|, + |A(0)|, = |A(^) - A(0)|, + |A(0)|, < e\z\^ + |A(0)|„ 
and so we conclude that 

(6.27) 1^1, <C|A(0)|^. 
One then observes from (16.251) and (I6.26P that for t > 0, 

(6.28) A(0)(t) = e^^Ui^(V'^(0)), 
where F : is given by 

F{v) := / e-^'^Us {Liv + N\v,e)) dr. 

J — oo 
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Note that F is evidently bounded and Lipshitz. We next claim that if ^/'(O) G M-e-, then 
-F('?/'^(0)) = 0. To see this, note that in this case the invariance property of the center manifold 
implies that G M.£ for t 7^ as well. Hence the unique fixed point of (16.261) must he z = Q 
and so in particular A(0)(t) = for t > 0, which establishes the claim. 

Finally, fixing any element ^/^^(O) G \ M.£ with sufficiently small iJ^-norm, one uses this 
last observation to obtain 

(6.29) ||F(^^(0))||^, = inf ||F(^^(0)) - F(^i)||^. < Cd^i(^^(0), A<,). 

tpi&Me 

The bound (I6.24p now follows by combining (I6.27p . (I6.28p and (I6.29p and using the bound 
(I6.23p . since z = ijj — ip- D 

Remark 6.6. Using the Implicit Function Theorem and the fact that is the identity on 
A^^, we find that is foliated on a small neighborhood of A^^ by transverse smooth manifolds 
P^^{w) through each w G Ai'^ , depending in a smooth fashion on the value of w. In particular, 
for (7 > sufficiently large and a > sufficiently small, the if ^-ball -8(0, a) is foliated by leaves 
P^^{w) for w G A^^ n 5(0, Ca), carried one to the other under the flow of the underlying 
ODE, uniquely specifled by the property that each solution initiating in P~^{wq) approaches 
the solution on the center manifold with initial data Wq at uniform exponential rate ~ 1 >> e. 
As a consequence, a stable manifold A/" C Ai^ fl 5(0, a) of an orbit or manifold of orbits 
within the center manifold Al^ extends to a stable manifold := VJyj(^j\fP^^{w) fl 5(0, a) 
in 5(0, a) of the same codimension in 5(0, a) as the codimension of M in Ai"^. That is, not 
only is asymptotic stability in 5(0, a) determined completely by asymptotic stability within 
the center manifold, but also conditional stability as measured by codimension of the stable 
manifold. 

We now apply the previous result on existence of a center manifold to assert the existence 
of bifurcating stationary and periodic states for equation fl6.12p . 

We begin with the case of stationary states bifurcating from the normal state. We refer 
to Section 4 for the deflnition (14.91) of the parameter xii which was found numerically to be 
real for I < Ic, positive for Ik < I < Ic and negative for < I < Ik where Ik ~ 10.93 and 
J, ^ 12.31. 



Stat 
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Proposition 6.7. For I fixed in the interval (0,1^), equation (16.121) exhibits a stable super- 
critical pitchfork bifurcation of stationary states {e^^°ilJe{-,£) '■ do G [0, 27r)} branching from 
the normal state for all sufficiently small and positive values of e. These equilibria satisfy the 
bound 

1 



(6.30) 



V-Xii' 

as predicted formally in Section 4 

For I fixed in the interval {Ik,Ic), the equation exhibits an unstable subcritical pitchfork 
bifurcation of stationary states e^^'^ipe{-,s), Oq G [0,2n), branching from the normal state for 
all sufficiently small and negative values of e. These equilibria satisfy the bound 

1 



(6.31) 



< C\e 



3/2 



as predicted formally in Section 4- 



Proof. Since we work here in the setting where / < Ic, the center subspace is spanned by 
the single eigenfunction Ui. Hence, we will express any point on the center manifold A4 as 
(<I>(/3, e),e) where /? G C corresponds to the coefficient of the point Pui on the center subspace. 
We begin with the case where / is fixed to lie in the interval (0, Ik)- Then for any small value 
/3o G C we let ip"^ denote the solution to (16.121) satisfying Dirichlet boundary conditions and 
initial condition tp'^{-,0) = $(/5o,£:). For all small, positive t, we know from Theorem 16.41 that 
{4'^{-,t),e) G Ai and so there exists a smooth function, which we denote by P"^ = P'^it), such 
that 



iq^ (6.32) ilj%x,t) = ^{f3'{t),e). 

Recalling the definition Pi := J^-^^ufdx we now apply the projection lie to every term in 
the equation 

satisfied by ip"^ and then integrate against Ui. We find that 

j^Ii,{rt)uidx = ^^ j Ii,{r)uidx=^^ j Yi,{^{P',e)uidx)=PiPl. 



stateode 



edefn 
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ac 
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We also have 

J Ilc{Li'ip'^)ui dx = Pi J Ui Liip^ dx = (3i j LiUiip'^dx = 

and J^^ Ilc{£'ip'^)ui dx = e[3i[3^ . Consequently, we obtain for [3^ the O.D.E. 

1 

PI = eP' + — Ii,{M[P'ui])ui dx + e(/5^ e) 
Pi J-i 

(6.33) = eP' + xnm'P' + e{P',e) 

where we recall the calculation of J^^ Ilc{-Af[P^Ui]^Ui dx carried out in Section 4, and we have 
introduced 

(6.34) e(/3^£) := l^,(AA[$(/3^ e)]) - Il,{j\f[P'ui]) . 

Pi 

Due to fl6.17p . we know that 

(6.35) e(e^^°/3", e) = e'^°e{P', e) for any e M 

and through elementary use of the triangle and Cauchy-Schwartz inequalities applied to the 
nonlinearity A/", along with (16.181) . we estimate 

|e(/3^ e) I < C m^P', e)] - ^[P^u^] ||^^ 

< C {mP^e)\\l, + mP',e)-P^u,\\^.) 

(6.36) = O {eiP'f + ip')') . 

Returning to ( I6.33p . consider first the case where P^{0) = Po & M.. We first claim that the 
function P'^it) must be real. To see this, we begin by noting that since Po is real, the quantity 
PoUi is PT-symmetric. Hence, in particular ip'^{-,ti) = ('?/'^(-, ti))^ as well for any fixed ti > 
since satisfies a PT-symmetric initial condition $(/3o'Ui,£:) that is . Now denote by ip^'^ the 
unique solution in to the equation 

(6.37) r'' = nP%h)uue,r'')- 

By (I6.32P we have ■ijj^{-,ti) = ^p^'^{■^0), and consequently, ^p^'^{■,0) = {'ijj^'^{-,0))^ . Evaluating 
fl6.37p at t = and applying the f operation to both sides, we then conclude that 

(6.38) i^''\-,0) = T{P%hrui,e,ij'''){;0) 



betaplus (6.39) 



closestart 
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as well. Applying the projection He to both (I6.37p evaluated at t = and fl6.38p . we see that 
indeed = as claimed. 

An easy application of the implicit function theorem then reveals the existence of a smooth 
curve of zeros e = e{P'^) to the equation 



(5^ 







such that e = —XniP^Y + ^{{P'^Y)- Hence, there exist smooth curves of equilibria P^{s) to 
f l6.33p for all small, positive e with 

1 



,1/2 



+ 0(6 



3/2N 



V-Xu 

Consequently, within the collection of points on the center manifold of the form with 
f3 real, the functions ipf := ^{/3^{e),e) represent a supercritical pitchfork bifurcation of 
equilibria from the normal state. The bound fl6.30p follows immediately from fl6.18p . In light 
of the rotation invariance of the problem, it immediately follows that there is in fact a circle 
of equilibria e^^°ipe, Oo ^ [0, 2tt) where we have written simply ipe for i^t- 

Regarding stability of these equilibrium, it is clear from (16.331) and the estimate (16.361) that 
given any initial data on the center manifold of the form $(/5oMi, s) with j3o real, positive and 
say bounded by C-y^, the solution to (I6.33P will converge to l3~^{s) and so the solution to 
(I6.12P will converge to ipe- Then since in light of (16.351) . (16.331) is clearly rotationally invariant, 
it follows that for complex initial data on the center manifold, i.e. initial data of the form 
^{PoUi,e) where Pq = |/?o| e*^" for some non-zero phase 6*0, necessarily the solution will converge 
to e*^°<l>(/5+(£:), e). Thus, one concludes that the circle of equilibrium states {e^^°ipe '■ 6o G 
[0,27r)} is asymptotically stable on the center manifold. 

Finally, suppose that we start with initial conditions for ( ]6.12p that are close to the circle 
of equilibria but that do not lie on the center manifold. That is, suppose we have 



(6.40) 



^^(0)-e'^»^e|Li <ri for some 6*0 G [0, 27r) 



but that '?/''^(0) ^ M-e- Without loss of generality, we ignore this rotation for the remainder 
of the argument and take Oq = 0. We will argue that for ri and e sufficiently small, again 
the trajectory ip'^it) is exponentially attracted to ipe- We take in particular, ri < where 



tritwo 



unstateode 
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Ci > 1 and r are the constants appearing in Theorem I6.5[ Since ||^e||//i ^ 
assert that 



V-xi 



41 

e we can 



(6.41) 



11^^(0)11^. <r/2 



by choosing e sufficiently small and appealing to (16.401) . Denoting by if)^ the trajectory on A^^, 
it then follows from (16.241) and (16.401) that for as long as il){t) obeys the bound ||^/'^(t) || < t, 
one has the estimate 



psmall (6.42) 



^-(t)-^^(t) 



< Cine-"* < -e""*. 
m 4 



Then the triangle inequality to implies that 



(6.43) 



< (1 + Ci)ri. 



Now il)'^{t) G M.£ is necessarily given by ^/'^(t) = $(/5'^(t), e) where 13^ is governed by (I6.33p . 
As we already noted, up to a rotation which we again ignore, the equilibrium value P^{s) that is 
stable under this flow, both in the sense of (exponential) asymptotic approach, /3^(t) 13^ {.^) 
and in the sense that (3'^{t) will stay close to this equilibrium for all time. Choosing ri still 
smaller if necessary, we may appeal to (I6.43P to conclude that |/5'^(0) — [3^{e)\ is small and 
then the Lipschitz property of the map $ allows us to assert that, for instance. 



tri3| (6.44) 



^^(t) - < r/4 for all t > 0. 



//I 



It follows from fICT]) . and (ICTj) that in fact Wip"" {t)\\ j^i < r for alH > and so 

is valid for all time. Combining (16.421) with the exponential approach of -ip^ to ipe along the 
center manifold, we obtain the asymptotic stability of all trajectories 'ip'^{t) satisfying (I6.40p . 

The case where Jfc < / < /c is handled similarly. Recall that in this case, the parameter 
Xii takes a positive value. Working then with e small and negative, and again starting with 
initial data of the form $(/?o,£^) with Pq real, we find that (16.331) is now replaced by 



(6.45) 



Another application of the implicit function theorem reveals that (16.451) possesses a pair of 
unstable equilibria /3^(e) for all small negative e-values with /^^(e) = ±^= kl^^^ + C^(kl^^^)- 

V All 
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An examination of fl6.45p shows that the corresponding equihbria ipf := e) bifurcat- 

ing subcritically from the normal state are unstable as well since even nearby PT-symmetric 
points on the center manifold, that is points of the form (3ui for (3 real and near I3^{e), flow 
away from them. Writing simply ipe for , the same is of course true for any of the equilibria 
on the circle {e*^o?/ie : 6q e [0,27r)}. 

□ 

We turn now to the case where the applied current I satisfies I > Ic and a bifurcation to a 
periodic state occurs. We recall that in this parameter regime the constants xii and Xi2 are 
not real. 



Hopf Proposition 6.8. Fix the applied current I in the interval (Jc, oo). Then provided that the 
constants x ct^d x given by and fl5.13p and fl5.14p respectively both have positive real parts, the 
equation fl6.12p exhibits a stable Hopf bifurcation to a periodic state ipp = iljp{x,t,e) branching 
from the normal state for all sufficiently small and positive values of e. This bifurcating 
solution, i/jp obeys the estimate 

(6.46) 11^, - e^^° {f3lp{t)v, + (3l^{t)v,) ||^, < Ce'/' 

for some 6q G [0,2tt), where the real functions 131 ^ and /?|p take the form 



stable2 



ampe 



albe 



Im Ai 



X cos 



VRex 
Imx 



0{e 



+ 0{e^'^) ) t 



3/2^ 



Rex 



sm 



ImAi + ^5 + 0(e^/2; 
Rex 



(6.47) 



Remark 6.9. Recalling the relationship between the functions Vi and V2 given by (16.90 . one 
checks that the coefficients /3f ^ and introduced above are related to the coefficients ai 
and a2 introduced in (15. 3p via the formulas 



ai{et) 



Aim Ai t 



''{(3lp{t) + t(3l^{t)), a2{et) 



-ilm Alt 



(Pipit) - tPl^^it)) 



Remark 6.10. We recall from the previous section that numerically, we indeed find that 
Re X > and Re x > for / > Ic- See Figure [61 
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Proof. Invoking Theorem 16.41 given any two complex numbers jS^ and of sufficiently small 
modulus, let ip"^ denote the solution to fl6.12p subject to Dirichlet boundary conditions and 
initial conditions given by <l>(/?5', /?2) Then we may describe via $ at all future times as 
= $(/3f (t), e) for complex-valued functions Pfit) and 
We then project fl6.12p onto the center subspace and use (16.101) to obtain 

— Im Ai/5|t;i + Im Ai/^^f 2 + ePfvi + e/3|t;2 

+ viAf [Pfvi + dx^ ^1 + ( / ^sA/" [Plvi + dx^ V2 
1 



+ yj ^vi i^MPl PI, e)] - Af[f3lv^ + (3lv2]) dxj v, 
+ { f V2 {M\^(5l /?2^ e)] - U[Plv^ + f5lv2]) d:^ V2. 




-1 

Integrating this equation first against vi and then against V2, we use the resulting two by two 
linear system in jSf ' and /5| ' to find 

Pt \ ^ I j',ViN'{PfVi+P'2V2) dx 

131 ) V-i ^2^f i/3fvi + PIV2) dx^ 

7_\ V, [UWl (31, e)] - ^^[Pfv, + PIV2]) dx' 
J_\ V2 {AfMPl (31, e)] - f/[(3lv, + (3lv2]) dx ^ 

Appealing to the center manifold estimate (16.181) and carrying out a lengthy calculation 
similar to that of Section [5], we finally arrive at a system of the form 

-[Rex{R'? + ilmx7^]/5f + [\mx{R'f - zRex7l/5| 

bi^ (6.48) +0 {e{R'f + {Wf) . 

In the system above we have introduced the notation 





(6.49) := J \(3l\' + \(3l\' and 7^ := z ((/3iT/3f - /^^^(/^l)*) = l/^i I l/^ll sin(^i - 



2)1 



where /3| = | /3| | e 
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We now apply the standard method for proving the existence of a periodic solution to (16 .481) 
via a Hopf bifurcation and to compute rigorously the amplitude and period of the oscillations. 
To this end, we consider first the case where the initial values Pi and /?2 real. Then the 
resulting PT-symmetric initial data $(/5°,/52,£) for fl6.12p will lead to the PT-symmetry of 
the solution at all future times. Consequently, for all t > 0, the projection of ip"^ onto the 
center subspace must take the form /3f + /5|(t)f2 where /3f and /3| are real. This leads to 
a significant simplification of fl6.48p in that 7"^ = 0. 

Converting to polar coordinates, R'^ and 6*^ := tan~^ {^2/^1)^ '^^ derive the system 



(6.50) {R^y = eR' -RexiR'f + 0{e{R')^ + {R'f) 



thetaeqnl (6.51) (^ )' = Im Ai - Im x{R')^ + O {e{R'f + {R'Y) 



Starting with any initial condition with small amplitude a := -R^(O), it is easy to argue that 
9'^ is a monotone increasing function of t, thus justifying the description of R'^ as i?^(a, 9^) via 
the scalar O.D.E. 

d p R p A/ 

^^^^ («'^^' = MAl^' - tat («')'+*(^- 

where 

g{e, R', 9') = O {e{R'f + {R'f) . 

A periodic orbit for ( ]6.48p corresponds to a value of the amplitude a such that i?'^(a, 27r) = a. 
Thus, we seek a fixed point of the Poincare return map 11 given by 
(6.53) 

n(e, a) := ae^'^'^^'^^' + e^WimA^ f ^-{e/imx^)e ( z5^i?-(a^ ef + gie, R^a, 9), 9)] d9. 

Jo V Im Ai J 

Using (I6.52P one readily checks that 



pomcare 



aeest | (6.54) R' (a, 9) = a + O^ea + a"^) and so R%a,9)^ = + 0{ea^ + a^). 



Therefore, expanding 11 for small e we see that 



n(£, a) =a + - — —a - - — — ^a"* + Oiae^ + ed^ + a^) 
Im Ai Im Ai 
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Via the implicit function theorem one derives a curve of zeros 6{a) for the expression 

Il{e, a) — a 
a 

Hence, we obtain a curve of fixed points of 11 with 

e{a) = {Rex)a^ + 0{a^), 

or inverting this relationship, 

(6.55) a = -^s'^' + 0{e). 

VRex 



Denoting the resulting periodic solution to (16.501) - fl6.5ip by (-Rp(t), dpit)) and letting {f3fp, /5|p) 
(-Rp cos(6'p), i?p sin(6'p)) denote the corresponding periodic solution to fl6.48p . the asymptotic 
estimates (I6.46P and (16.470 then readily follow from (I6.5ip . (16.540 and (I6.55p . 

We also note that the asymptotic stability of this periodic orbit, and hence the asymptotic 
stability of the periodic solutions to (]6.12p among nearby PT-symmetric competitors on the 
center manifold, is a consequence of the condition 

•^n . . . . 47rRe X 2 ^ ^ 

— (e(a), a) ~ 1 - - — —a < 1, 
oa Im Ai 

cf. [6], Thm. 12.13. In light of the rotational invariance of the system (16.480 . this means not 
only that real initial data (/3i(0), /3|(0)) lying sufficiently close to the periodic orbit will be 
asymptotically drawn by the flow into the orbit but also that any rotation, say e*^°(/3f (0), /?f(0)) 
of such initial data will be drawn to the corresponding rotation of this orbit as well. 

We next wish to argue that these periodic orbits are stable within the class of all flows 
starting nearby on the center manifold. For this purpose we return to (16.481) : that is, we 
consider the situation where the initial data for jSf and /5| are not necessarily real. Note that 
this corresponds to initial data that is then not assumed to be PT-symmetric. We claim that 
if we start with complex initial data (0), /3|(0)) sufficiently close to any of the orbits given 
by (16.460 - 06.470 . then again the flow will draw the resulting solution into one of these periodic 
orbits. 

To establish this claim it is useful to derive a system of O.D.E.'s for the quantities (-R^)^ 
and 7"^. Differentiating the deflning formulas in (16.490 and using the system satisfled by /5f 
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and /3| given in (I6.48p . a routine calculation leads us to: 

(6^56) M' /2[(.-Rei -Re,W]\ ^ ^(^j^.,. , 

\ 'l' / \ 2\e-Re(x + xm?]Y ) 

Note that /3| are parallel if and only if 7^ = 0. Recalling that PT-symmetric solutions cor- 
respond to /5J real, we thus see that 7^ = corresponds to the invariant three-dimensional 
manifold of PT-symmetric solutions and their complex rotations that we have already consid- 
ered, with the single remaining dimension parametrized conveniently by 7^. 

If we momentarily ignore the error term, then we can readily identify an equilibrium point 
for fl6.56l) at (e/Rex, 0). Linearization about this critical point immediately yields linear 
stability and in fact a standard phase plane analysis yields nonlinear asymptotic stability. To 
treat the full system fl6.56p (that is, including the error term), we must linearize instead about 
the nearby periodic solution ((i?p)^,0), corresponding to the (real) periodic orbit (/3i_p, /?2,p) 
constructed above, yielding 

{Y)' = 2 [5 - Re (x + X){RI? + 0{e{Rlf + (i?;)^)] 7^ 

Since the error term is clearly lower order, the condition {RpY ~ e/Rex again implies ex- 
ponential linearized stability in the Y direction. As exponential orbital stability has already 
been verified in the remaining directions, we may conclude exponential linearized and non- 
linear orbital stability of the family of periodic orbits within the full four-dimensional center 
manifold and not only within the invariant three-dimensional 7 = manifold. 

It remains to argue that the periodic orbits attract nearby initial data for fl6.12p off of the 
center manifold. The argument for this fact follows exactly as did the corresponding point in 
the proof of Proposition 16.71 regarding stability of stationary solutions through an appeal to 
Theorem 16.51 □ 

Proposition 16.81 describes the bifurcating stable periodic solutions of main physical interest. 
With a little further effort, we may also rigorously obtain the existence of the two unstable 
periodic branches derived formally at the end of Section [5] and obtain essentially a complete 
description of the full four-dimensional bifurcation. 
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Hopf2 Proposition 6.11. Under the assumptions of Proposition {67^ the equation (16.121) also ex- 
hibits an unstable bifurcation for small e > to periodic states ipp^± = 4'p,±{x,t,e) of form 



uiistable2 



(6.57) 



where, recalling that Xu 



1 



ix + X) so that — Rexii > 0, 



unstable3 1 (6.58) /3i'^(t) 



+ C(r 



3/2\\ ±j({ImAi+eij^)i+0{£3/2)) , ne,± ,-ne,± 

' )]e ^ Hcxii' V ana P2 = ±ipi . 



Furthermore, there exists a positive value a independent of e such that for all small e > 0, 
these two states, along with the stable periodic state ipp constructed in Theorem \6.^ and the 
normal state, represent the only persistent states in B{0,a) := {tp : < a}. Indeed, the 

phase portrait within B{0, a) consists of two repelling codimension two C°° stable manifolds 
Mp^± of the unstable periodic solutions ipp^±, a repelling codimension four C°° stable manifold 
A/q of the normal state ip = 0, and an attracting codimension one C°° stable manifold MpT 
of the invariant manifold Mpx of PT-symmetric solutions and their rotations lying within 
Ai"^ n B{0, Ca) , cf. Remark \6.(k The latter contains both the unstable normal equilibrium and 
the stable periodic solutions, with all other solutions flowing from repelling submanifolds Mp^± 
and A/q to the attracting submanifold Mpt and thereafter to the stable periodic solutions. 

In particular, all solutions originating in B{0,a) converge either to the stable periodic, an 
unstable periodic, or the normal state. For generic initial data tpQ G -8(0, a), the corresponding 
solution converges to the stable periodic state, the only exceptional data lying on the codimen- 
sion two and four submanifolds M'p^± and Mq, respectively. 



Remark 6.12. The description (I6.57p - (l6.58p of the unstable periodic solutions given in 13- 
coordinates may be recognized as profile (I5.16p - (l5.17p given in the a-coordinates of Section 

El 



Proof. Our approach here is generally to establish the claims of the theorem first for the 
system of differential equations (16.481) or (16.561) ignoring higher order error terms and then to 
broaden the claims to the full equations including the perturbations. Hence, we first note that 
in light of the smoothness of these higher order terms and the polynomial error bounds, the 
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perturbations are seen to be small in the topology in a small neighborhood of the origin, 
i.e. for j3l, /5| small or equivalently, for small. 

We begin our analysis by turning to the system (16.561) . For convenience, let us denote 
the quantity (-R^)^ by and so view this system as one for 7^ given by (I6.49P and A'^ : = 
+ l/?lr- We list some simple observations: 

• Elementary phase plane analysis of this system reveals that there exists a small positive 
number r, independent of e, such that the triangle 

T:={(7^A^) : 0< |7"| < A' < r} 

is invariant for the flow (16.561) for all e sufficiently small. That the flow cannot exit the two 
bottom sides of this triangle is trivial in light of the algebra 

h'\<2\(3l\\f3l\ < \(3t\' + \(3l\' = A^. 

That the flow cannot exit the top where A"^ = r follows immediately from the fact that 

(6.59) {A')' < 2(^6 -Re x{A') + Ce{A') + C{A'Y^ A' < 

provided that e << A^ and = r is small. Here C is a positive constant coming from the 
error bounds in (I6.56p . 

• A further consequence of (I6.59P for the full system (I6.56P and its "parent" system (16.480 
is that if one starts with data lying in the triangle T, then the corresponding orbit will 
exponentially approach a sub-triangle in which A"^ = O^e). The rate of approach may be 
obtained by direct computation from the differential inequality (^4"^)' < — Ci(y4^)^, Ci > 0. 

• If one ignores the error terms in this system, then it is easy to check that there are exactly 
four critical points located at 

(7^ An = (0, -^), (0, 0) and ' (±1, 1) 
Kex ^l^^-exiil 

and the Jacobian at each of these points is nonzero (see the next item below). 

• Linearizing (16.561) about these four critical points, one finds: (i) The critical point near 
(0, ^f^) is stable, with two negative eigenvalues, as has been noted already in the proof 
of Theorem 16.81 This point corresponds to the stable periodic orbit arising from a Hopf 
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bifurcation, (ii) The normal state (0, 0) is unstable with a multiplicity two positive eigenvalue, 
(iii) 2|Re^xiil both saddles with one positive and one negative eigenvalue. The 

negative eigenvalue corresponds to the eigenvector (±1,1). 

• Not only do (±1, 1) correspond to the stable eigendirections of the linearization about 
2|Re^xiil respectively, but in fact, the bottom sides of the triangle T, i.e. Y = 
exactly the stable manifolds of these unstable critical points for the system (16.561) ignoring 
error terms. 

• Introducing s as the ratio 7^/^"^, we use (I6.56P to calculate that 
(6.60) s' = -2Re xA's{l - s^) + 0{eA' + {A'f) . 

As a consequence, if one ignores error terms, then one sees the flow repels away from the 
triangle's lower sides s = ±1 towards the manifold of PT-symmetric solutions and their 
rotations, s = 0, i.e. towards 7^ = 0. We should note that since the flow (16.121) preserves PT- 
symmetry, it follows that s = must in fact be a critical point of (I6.60p even with inclusion 
of error terms. 

Having noted these simple properties of (16.561) . we next establish the existence of the unsta- 
ble periodic orbits for (16.481) and hence for (16.121) . As we shall see, these correspond precisely 
to the two unstable critical points near 2pf^^(±l,l) found above in the 7^^^ plane. Let 
us first note that back in Section [5], we identified exact periodic solutions on the invariant 
manifolds ai = and a2 = for the equations (I5.6l) - (l5.7p obtained by neglecting higher-order 
error terms. In light of the relations ( 16. 9p this corresponds to exact solutions (again neglect- 
ing error terms) to (16.481) found along the manifold /3| = ±i/3^. This, in turn, corresponds to 
solutions of (I6.56P for which 7^ = ±A^. 

The explicit formulas for these two unstable periodic solutions can be obtained by directly 
solving (16.481) under the constraint /3| = ±2/5^ and are given in (I6.58p . Note in particular 
that the two equations in (I6.48P for (Pf)' and (/3|)' reduce to the same equation under this 
constraint. Their saddle- type instability follows from the observations listed above and in 
particular from (I6.60p . 

In order to argue that these saddle-type periodic solutions persist with the inclusion of error 
terms in (16.480 . we define cjf := (/3f + z/?|) and tuf := (/?f — i/^f) so that cuf and (-jf are closely 
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related to the aj variables of Section 5 via the formulas cu^ = 
We will describe how to find a periodic solution to the full system fl6.48p nearby the solution 
corresponding to jSf = (i.e. c^l = 0) for the truncated system. The same analysis can be 
repeated to find the other unstable solution to the full system as well. 

To this end, it is perhaps best to introduce a Poincare map as follows: Consider (16.481) 
without error terms and start with initial data c^KO), t^|(0) such that Imco'f(O) = while 
Re (0) is near ^/e/ v^— Rexii and |ci;|(0) | is small. We then measure the values of Re cuf , Re cul 
and Imcuf at the next time, say T = T(c<j^(0), c<j|(0), e), at which the trajectory (cuf (t), co'l(t)) 
crosses the 3-plane {(ci;i,co'2) gCxC: lmco'i = 0}. Let us denote this map 



(6.61) 



(Reu;f(0),Reo^|(0),ImL^^(0)) ^ (Reu;f (T), Reo^^(r), Imu;^(T)) 

Note that (^=^=,0,0) is a fixed point of this map corresponding to the 



by F= : R3 - 

exact unstable periodic solution of the truncated system satisfying = -i/^f. 

By converting fl6.48p into a system for Uj, j = 1,2 one can carry out a tedious but straight- 
forward calculation to determine that when the Jacobian matrix of is evaluated at the 
fixed point ( ^_^^^^ , 0, 0), it takes the form 



(6.62) 



Here 



and 



DF' 











e'* cos 6^ — e" sin&^ 
e""" sin ¥ e"'" cos ¥ 



Rex 
-Rex 



11 



27r 



Im Ai 



Imx 



12 



Rexii 



ImAi + 



Rcxii 



which we recognize as the period of the unstable periodic solution with error terms ignored. 

Note that there is one eigenvalue of DF"^ of modulus less than one, namely /i^, giving one 
stable direction to the map, while the other two unstable eigenvalues e"''^*^^ have modulus 
greater than one since > 0. Thus, one again sees the saddle- type instability of this periodic 
solution from the perspective of this Poincare map. 
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To argue that this picture persists under perturbation (i.e. under inclusion of the error 
terms in fl6.48p ). we observe that the three eigenvalues of D{F^ — I) take the form 



Once we incorporate the error terms from (16.481) in our analysis, we use the fact that within the 
0{^/£) ball where we are working, these (9(£|/5|^ + |/5|^) error terms are necessarily 0{e^^'^) with 
a contribution to the entries of the original Jacobian matrix DF'^ of at most (9(e|/?p + = 
0{e^). Hence, again subtracting the identity to form the displacement map of the full system 
including error terms, (I6.63P guarantees that the differential of this map is still nonsingular. 
Necessarily it then must still map a neighborhood of the original critical point onto the origin, 
that is, there must still exist a fixed point of the map corresponding to F perturbed by higher 
order terms. What is more, the Jacobian matrix of this perturbed map evaluated at the fixed 
point must still possess one eigenvalue of modulus less than one and two of modulus greater 
than one so the saddle-type instability persists under perturbation for these periodic orbits. 

It remains to establish the asserted global behavior of solutions originating in the H^- 
ball B{0,a). We begin by describing the behavior in A^^ fl 5(0, Ca), where C > is the 
geometric constant of Remark 16. 6[ to be used later. As observed previously, the manifold 
A/pT C Ai^ of PT-symmetric solutions and their rotations remains invariant for the full as 
well as the unperturbed equations, and is locally attracting by (16.601) . Likewise, the normal 
state {PI = /^l = 0) is an equilibrium of the full system that is locally repelling. Since the 
unstable periodic solutions correspond to hyperbolic, saddle-type equilibria of the associated 
return map, with a single stable eigenvalue and a pair of unstable eigenvalues, we find that 
they possess, for fixed e, a C°° stable manifold of dimension one of the return map. Under 
time-evolution, this induces stable manifolds A/^,± of dimension two within the center manifold, 
as claimed. 

The structure and attracting or repelling properties of A/^,± are difficult to determine outside 
an e- neighborhood of the periodic solution, due to the fast, order one, angular flow relative to 
the order e flow measured by the return map. Thus, it would appear that a determination of 
the global structure in Ai"^ fl B{0, Co) would require more complicated averaging arguments 



smalleig (6.63) 



Im A 
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outside the scope of the present analysis. However, we may finesse this point using more 
elementary tools together with the special structure of our equations. 

We first observe by (16.601) that for C > sufficiently large, the neighborhoods K,± of the 
triangle T given by 

Kd^ (6.64) /C± := {(7^ A')eT: s := Y/A' satisfies \s±l\<C{e + A')} 

are invariant under backward fiow of the ODE on the center manifold. Hence, by estimate 
(I6.57p -( j6758l) showing that — 1| < Ce along the unstable periodics, the stable manifolds Ap,± 
of the unstable periodics are confined to IC±. In particular, they lie no more than C{e + A'^) 
from the sides of the triangle corresponding in j^-A"^ coordinates to the stable manifolds of 
the unstable periodic solutions of the unperturbed equations. Moreover, from (16.561) . we see 
that for (7^, A'^) G IC± one has 

{A')' = 2{e - Re (x + x)A')A' + 0{s + A'){A'Y. 

Thus, for (7^, A^) G /C± outside an ^^-neighborhood of the critical point p^^-^4-j-^(±l, 1) of the 
unperturbed '-y^-A'^ system, is strictly decreasing for A^ > (j^+x) ^^'^ strictly decreasing 
for A^ < (ji+x) ' ^^^'^^ these regions, A/^,± are graphs over their unperturbed counterparts 
and 0{e + A"^) close to them. 

It remains to treat the excluded neighborhood in 7^-y4^ coordinates of the critical point 
^j^4q-^(±l, 1) of the unperturbed j'^-A'^ system. By a straightforward computation, we find 
that this neighborhood is contained in the image of an e-neighborhood in /^^-coordinates of 
the orbit of the order e^^"^ amplitude unstable periodic solution (16.581) of the unperturbed 
/3^-equations. 

Indeed, a brief examination in /^'^-coordinates shows that the ^^-A^ estimate is overly con- 
servative (a result of the singularity in certain directions of the coordinate change between 
the two coordinate systems). Changing to the more convenient cj-coordinates, for which 
1^1 1 = ^J^^xii ' '^2 = corresponds to the unstable periodic orbit, and tuf = its stable 



Nbd 



Thetabd 
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manifold, we find compute that 

\ul\t = {e + (Rexiilcu^P + Re Xi2\^l?)\ujl\ + 0{e\uj'\'' + 

\ujl\t = {e + (Rexii|cu||^ + Re xi2\i^{\^)Wl\ + 0{e\uj'\^ + 

(cf. f l5.10l) -( l5.11l) ). Using this system, we can check that as long as the deviations uj\ and a;| 
from the unstable periodic solution to the unperturbed system satisfy 

£3/2 « p^l « ^1/2 and e""^ «\Col\«e'l\ 



the stable manifolds Ap,± are graphs over their unperturbed counterparts 



in this case) 



and all other solutions are repelled toward MpT since ^ | will be negative in this regime 
while ^ |(2;|| will be positive. 

On the other hand, by our previous estimates on the orders of perturbation terms and their 
Jacobians in the return map F"^ given by (]6.6ip . the Taylor expansion of this map at the 
unstable fixed point corresponding to the unstable periodic, denoted here by p^, is 

F^z) ■= F%pl + z)-pl = dF%pl)z + N'{z, z) + Q%z). 

Here the quadratic order Taylor remainder term N"^ is order 

(6.65) N' = Oi\pl\\z\'), N',=Oi\z\' + \pl\\z\) 
and the term which incorporates the perturbation term is order 

(6.66) = Oielff + \pl\'), Ql = Oie\ff + 

in an e-neighborhood of pi- 

Reviewing the standard invariant manifold constructions by fixed point/contraction argu- 
ments (see, e.g., [3|, [6]), we find that they yield existence and closeness in angle of these 
manifolds to corresponding invariant subspaces of dF^ on a ball about pi for which the Lip- 
shitz norm of the (total) nonlinear term A^^ + 6^ is sufficiently small compared to the spectral 
gap of dF'^{0) = dF'^ipl) so long as the norm of log dF'^ = logdF^ is no larger than some 
specified multiple of the spectral gap, as it is here, cf. f l6.62p . Here the spectral gap is defined 
as the minimum distance between one and the modulus of eigenvalues that are not modulus 
one. By (16.631) . the spectral gap is greater than rje for some positive r], while by fl6.65l) - fl6.66l) 



54 J. RUBINSTEIN, P. STERNBERG AND K. ZUMBRUN 

and the fact that \pl\ ~ e^/^, we have |A^|| + |9^| < C(|zp + e^^'^\z\ + e"^) = o{e) as desired 
for \z\ « e^l"^ . Thus, we obtain a detailed "microscopic" description of the behavior of the 
return map on a ball about of radius r/e:^/^, for 77 > and sufficiently small. This trans- 
lates to a detailed description of A/^,± and asymptotic behavior on an ?7e^/^-neighborhood of 
the unstable periodic orbit, far more than what was needed (the excluded e^/^-neighborhood 
in u;^ coordinates or for that matter the excluded neighborhood in 7^, coordinates) to 
complete the argument. 

Indeed, though we do not need it, the faster decay in e of perturbation terms O, yields 
uniform convergence of perturbed to unperturbed flow and stable manifold in "blown-up" 
coordinates := cj/e^/^, by the same argument. 

Thus, solutions originating outside /C± remain outside, with s by f l6.60p strictly decreasing at 
an exponential rate. Moreover, if ever solutions leave /C±, they are attracted at an exponential 
rate to the PT-symmetric manifold Mpt corresponding to s = 0. On the other hand, solutions 
remaining in 1C± sufficiently long must eventually enter the small neighborhood 5^ fl /C-t, 
after which solutions not on A/^,± must (by our microscopic description carried out in lo- 
coordinates) eventually leave. Piecing together this information, we find that all solutions 
originating in A^^ fl 5(0, Co) and not on the stable manifolds Ap,± of the unstable periodics 
time-asymptotically approach the attracting PT-symmetric manifold Mpt^ as claimed. This 
completes the description of asymptotic behavior within the center manifold. 

By Remark 16.61 the stable manifolds A/^,± within the center manifold, of codimension 
two in hA'^ n -8(0, a), extend to C°° stable manifolds 

iVp,± := U^eAr,,±P,"'(w) n 5(0, Ca) 

of codimension two in 5(0, Ca). Consequently, for fixed e > 0, (Ap,+ U Ap _) fl 5(0, a) is 
exactly the set of data in 5(0, a) whose solutions converge asymptotically to an unstable 
branch. Likewise, there is a stable manifold A/q := Pg^^(O) of the normal equilibrium, of 
codimension four, containing all solutions originating in 5(0, a) and converging to the normal 
state. Finally, Mpt has a C°° stable manifold Apt := U^gA/'pr-Pe H'"^) of codimension one. By 
Proposition 16. 5[ together with our description of asymptotic behavior on the center manifold, 
we find that all solutions originating in 5(0, a) outside the sets Ap,± and Ao are attracted 
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to NpT (hence to the larger manifold A/pT of solutions converging to Npr) and ultimately to 
the stable periodic solutions or their rotations. This completes the description of asymptotic 
behavior and the proof. □ 

Remark 6.13. An implication of Proposition 16.111 is that if one chooses complex initial 
data (/5^(0), /5|(0)) away from the special two-dimensional manifolds containing the unsta- 
ble branch, then there is a kind of phase-locking phenomenon whereby the flow (16.481) pushes 
the solution towards a pair (/5i(i), /5|(t)) which is a complex rotation of a real pair, that is, 
7*^ tends to zero and the corresponding solution to (16.121) approaches a PT-symmetric profile. 
On the other hand, it is not true that all solutions exhibit this phase-locking, as evidenced by 
the existence of the unstable periodic solutions. 

Remark 6.14. As noted parenthetically in the proof, a side-consequence of our analysis is to 
rigorously validate the formal asymptotics of Section [5] by verifying convergence of perturbed 
to unperturbed phase portrait on a ball of size e^/^ about the normal state, in "blown-up" 
coordinates Cj := uj/e^/'^ (equivalently, (3 := P/e^/"^) equivalent to the a-coordinates of the 
earlier section. 

6.3. Bifurcation off of higher eigenvalues. We conclude this section by remarking on the 
bifurcation situation off of higher eigenvalues. Recalling the notation fl6.3l) . we fix any integer 
£ > 1, and set the parameter F in (11.41) equal to Re X2i-i+e. Then f ll.4l) - fll.5p can be rewritten 

as 

(6.67) iJt = U[i,]+M{iP,e) 

All of the analysis of the preceding subsection, including the center manifold construction, 
applies to produce the existence of stationary and periodic solutions to this system for I < h 
and I > Ii respectively, cf. (16. 2p . The difference here is that for £ > 1, there will always 
be a nonempty unstable subspace corresponding to \j with 1 < j < 2i — 1. Hence, these 
bifurcating solutions will always be unstable. As far as the center manifold construction is 
concerned, the primary change is that the fixed point argument must now be applied to the 
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integral equation 

y{x,t)=V{u,e,y) := e^^^u + f e'^'^'-^^Il.U' {y{x,T),e) dr 

Jo 

/t poo 
e^^(*-^)n,Ar^(y(x,r),£) dr- / e^^(*~^)n„ A^-^ r), e) dr, 
oo Jt 

rather than fl6.19p . Here Ilu denotes the projection onto the unstable subspace and through 
its realization clS cl contour integral 



Jt-u 



"bd2] (6.68) e 

(where r„ is any bounded contour enclosing the finite number of eigenvalues with positive real 
part) one establishes the necessary bound 



|e*^^n„||^,^^i < Ce'^^'-'^' for alU < 



to augment flO^jl - flQajl . 



7. Phase slip centers 



The physics literature associates periodic solutions with the existence of phase slip centers 
(PSCs), that is, zeros of the order parameter. Indeed, an immediate conclusion of the formal 
calculation of Section [5|, rigorously confirmed in Section 6, is that the leading order term in 
the expansion for the periodic solution ipp to fl6.12p along the bifurcation branch in the regime 
I > Ic has a periodic array of zeros at x = 0. More precisely, referring back to ( ]6.46p -( !6.47p . 
we see that the approximate solution Pi{t)vi{x) + P2{t)v2{x), when evaluated at a; = 0, has 
zeros whenever the expression 



cos[ti;^t] := cos 



Rex 



vanishes, since ^2(0) = i(ui(0) — ^2(0)) = 0. In other words, there are PSCs located period- 
ically at {x,t) = (0,T|), /c = 0, 1,2, . . . where ~ ^(7r/2 + kn). Since the actual solution 
ipp is uniformly close to (3iVi + /52f2, it follows that the two functions have the same Brouwer 
degree in a neighborhood of the points {(0,T^)}. Hence, we will rigorously conclude that 
ipp possesses an array of zeros close to the points {(0,T^)} once we check that the degree of 
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PiVi + l32V2 is nonzero at these points. But a direct calculation of the Jacobian determinant 
at these points indeed reveals that 

Jac ipivi + P2V2) (0,T|) Reu[{0). 

A straight-forward numerical calculation of the first eigenfunction for the operator M (cf. 
(11. 9p ). for which there are rigorous error bounds, shows that ReM'i(O) 7^ 0, with values ranging 
monotonically from about —0.2234 for / = 12.5 down to about —0.3578 for / = 20. Hence, 
we obtain a rigorous confirmation of a periodic array of PSC's for the periodic solution. 
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